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PREFACE 
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This  monograph  is  basically  the  outline  of  lectures  given  at  the 
WECOH  Research  Seminar  on  Applied  Mathematics  and  Continuum 
, Mechanics  held  at  Watervliet  Arsenal,  New  York,  in  June  1973« 

h 

I The  purpose  of  the  lectures  was  to  familiarize  theoretically  ^ 

oriented  engineers  with  some  mathematical  techniques  which  are 
useful  in  continuum  mechanics.  Abstract  differentiation  and 
integration  is  shown  to  lead  directly  to  the  formulation  of 
known  or  new  variational  principles.  Some  applications  are 
also  given  to  specific  types  of  ordinary  differential  equations 
which  occur  in  theoretical  mechanics.  However,  the  formulation 
of  variational  principles  for  boundary  value  problems  of 
ordinary  differential  equations  can  obviously  be  extended  to 
a much  broader  class  of  problems. 
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Chapter  1 


Abstract  Differentiation  and  Integration. 

1.0  Introductory  conunents . 

In  this  chapter  we  shall  consider  primarily  problems  posed  in  Hil- 
bert space  setting.  Most  statements  offered  here  are  easily  trans- 
lated into  Banach  space  terminology,  where  we  have  to  keep  careful 
count  of  which  elements  belong  to  the  space  B and  which  belong  to 
its  dual  3*,  or  even  to  B**.  In  the  Hilbert  space  setting  we  can 
afford  to  be  sloppy,  and  the  arguments  are  frequently  sinplified. 

1 . 1 The  concept  of  a derivative . 

To  offer  generalizations  of  the  concept  of  a derivative  and  of  the 
usual  necessary  condition  for  the  extremum  of  a function  we  need 
to  have  a look  at  the  concept  o^  a derivative  of  a function  of  a 
single  (real)  variable,  and  of  two  variables.  Differentiability  of 
a function  f whose  domain  is  D(f),  f(x):  D(f)  c R R,  can  be  stated 
as  follows;  given  Xq  € D(f)  and  given  h,  we  can  express  the  dif- 


ference f(xQ  + h)  - f(Xp)  = Aj^f 


in  the  form 


(1) 


A jjf  = K(  XQ,*h)  *h  t »£.(  XQ,-h)  , 


where  K is  some  constant  (depending  on  Xq  and  on  h)  and  K(xQ,h) 
has  the  property 


lim. 

h-^0 


IhI 


= 0. 


Moreover  (regarding  K as  a function  of  h)  liir  K(x_;h)  exists. 

h->0 

We  could  of  course  generalize  this  by  considering  lim  , or  lim 

h-^0^  h,0_ 
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% , 


V.? 

¥ 


or  lim  sup  K(x-;h),  etc.  obtaining  the  Dini ' s right,  left,  upper 
h-*-0  ^ 


and  lower  derivatives. 

We  shall  deliberately  stay  with  only  the  siirplest  concepts. 
The  relation  (1)  can  be  rewritten  as 


= <K(XQ;h),h>  +t(Xp,h). 


where  in  one  dimension  the  inner  product  <,  > is  a simple  multipli- 
cation. 


Af 


At  this  point  we  have  deliberately  avoided  the  form  , 


since  generalizations  of  this  form  are  possible  only  in  spaces  where 
products  of  this  type  are  defined.  In  general  h will  be  a vector 


in  an  infinite  dimensional  vector  space,  f may  be  some  mapping  and 
Af 


the  product  — makes  no  sense  at  all.  However,  the  analog  of 
formula  (1)  has  an  easy  interpretation.  It  is  clear  that  for  the 


infinite  dimensional  case  one  dimensional  concepts  are  inadequate. 


and  we  shall  consider  the  concept  of  a derivative  in  R , and  attempt 

.2 


to  generalize  the  basic  notions  from  R to  an  arbitrary  Hilbert  space. 
We  shall  consider  two  basic  pointwise  maps 


f : 


2 1 

<|):  R -»■  R . 


f will  be  called  an  operator,  a map,  or  a transformation,  (|)  will 
be  called  a functional.  Clearly  a functional  is  a map  R^  ■*  R^((^ 


However,  the  one  dimensional  range  of  this  map  sinplifies  many 
arguments,  and  such  maps  should  be  considered  separately. 


! I 


'•an 


A 


3 


2 

The  map  f may  have  the  following  types  of  derivatives  in  R 

i)  partial  derivatives 

ii)  total  derivative  - or  the  Jacobian  derivative, 

iii)  directional  derivative  (in  the  direction  of  some  two  dimension- 

2 

al  vector  in  R ) . 

For  the  functional  $ we  introduce  the  concept  of  a gradient  of 

All  these  concepts  have  appropriate  generalizations  in  Hilbert 

space.  We  can  also  introduce  the  appropriate  concepts  of  integration 
. 2 

in  R . Given  any  simple  arc  T,  whose  endpoints  are  Pj^,  P2»  we  can 
introduce  the  Riemann  integral 

I <f(x,y),ds>  If  f(x,y);  R'^  -*■  R^  has  the  property  that  for 


2 

any  Pj^,  P2  ^ ^ ^ such  integral  is  independent  of  the  path, we  in- 
troduce the  idea  of  a potential  functional. 


(x,y) 


<MxQ,.y_)  = 0,  $(x,y)  = / <f(x,y)  . ds> 

( Xq  , y 0 ) 
r 

We  note  that  the  independence  of  path  condition  requires  no  proper- 
ties apart  from  Riemann  integrability . 

The  more  commonly  used  conditions  curl  f(x,y)  = 0 in  R^  or 


<Sfj^(x,y)  6f2(x,y) 


IT 


'S  X 


in  R 


f = (fj,  f2> 


assume  dif ferentiabili ty  of  f,  and  assume  siitple  connectivity  of 
. Each  of  these  criteria  will  be  shown  to  have  inportant  gener- 
alizations in  Hilbert  space. 

1 . 2 Definitions  of  derivatives  of  maps  in  Hilbert  spaces . 

By  a map  in  Hilbert  space  we  mean  a map  (transformation,  operator...' 


4 


f : n = H2 


where  Ji  is  some  subset  of  a Hilbert  space  and  the  range  of  f 
is  a subset  of  a Hilbert  space  H2 . 

The  simplest  case  when  H^^  = •fc  is  easily  disposed,  since  division 
by  a scalar  is  defined  in  H2.  We  say  that  f:  R -►  H2  is  differen- 

tiable at  the  point  tQ  € R if 


lim 

t-^t„ 


f(t)  - f(tp) 
t - t« 


exists . 


Unless  otherwise  qualified  lim.  will  stand  for  limit  in  the  norm 
(strong  limit).  That  is  we  postulate  the  existence  of  a vector 
f'(tQ)  € H2f  such  that 


lim 

t-^t. 


f(t)  - f(tQ) 


t - t. 


- f'(to) 


= 0, 


(m) 


We  can  rewrite  this  in  the  form  (1.1)  as  an  equation  in  H2 


f(t)  - f(tQ)  = (t  - tg)  f'Ctg)  t'KtjtQ)  (1.2) 

= 0,  (1.3) 

where  we  deliberately  avoid  complications  at  this  point  by  insisting 
on  this  limit  being  defined  instead  of  possible  alternate  conditions 
(one  sided  limit,  lim  sup,...) 

If  we  replace  convergence  in  the  norm  by  wealc  convengence,  we 
have  the  equivalent  definition  of  a weak  derivative: 


where 


lim 

t->t„ 


t - t„ 
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f;  R ■+•  H2  has  a weak  derivative  at  the  point  tg  € r,  if  there  exist 
a vector  f'(tQ)  ^ H2  such  that 

f(t)  - f(t.) 

f'(t-)  = weak  limit  ( ) , (1.4) 

t * t„  t - ‘0 

with  the  formula  (1.2)  still  valid,  and  condition  (1.3)  replaced  by 

<)^(t;t„)  w 

( -piz ) 0,  (1.5) 

^ -0 

where  0 is  the  zero  vector  in  These  definitions  don't  make  much 

sense  if  f is  an  operator  f:  ^ ^2  (generally  is  infinite 

dimensional) . 

The  definition  of  directional  derivative  due  to  Gateaux  general- 

2 

xzes  the  concept  of  a directional  derivative  in  R . We  consider 
the  map  f:  n = ->  H2.  Take  Xq  - an  interior  point  of  n.  Pick  a 

fixed  vector  h <;  H^.  For  fixed  Xq  , h ^ consider  the  vector 
(f(xQ  + th)  - f(xQ))  ^ H2  where  the  constant  t is  picked  suffic- 
iently small  in  absolute  value  to  make  sure  that  x^  + th  € which 
is  possible  since  Xq  was  an  interior  point  of  For  fixed  Xq  and 
h the  difference  <^f(t)  = (^(Xq  + th)  - £(Xq))  is  a function  of 

the  real  variable  t only.  Hence  we  are  back  in  the  previously 

discussed  case  5f(t);  R H2 , and  we  can  define  the  directional 

derivative  of  f in  the  direction  of  h,  computed  at  Xq , if  there 

exists  a vector  f(xQ;h)  € e^  such  that  for  all  sufficiently  small 
values  of  t,  the  formula  holds: 

<Jf(t;XQ,h)  = f(xQ+th)  - f(Xg)  = t*f'(xQ;h)  +'t(t;XQ,h)  (1.6) 


where  lim 
t->0 


0. 


, -S€e  [to]  ) . 


(1.7) 
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»t(  t;  Xgh) 


X w 

(or  -£\l(t;XQh)  0 as  t ^ 0 xn  case  of  a weak  derivative.) 

We  should  point  out  that  the  directional  derivative  f'(xQ;h)  does 
not  have  to  be  linear  in  h.  In  general  f ' ( Xq  ;hj^+h2)  f'(xQ;hj^)  + 
f (xQ;h2).  It  is  however  homogeneous  of  degree  one  in  h,  as  can 
be  easily  checked  from  the  definition,  and  f'(xQ,'ch)  = clUx^rh). 

1 . 3 The  Gateaux  derivative . 

We  now  assume  that  f'(xQ;h)  exists  and  is  linear  in  h.  Then  f'(xQ;h) 
is  called  the  Gateaux  derivative  of  f in  the  direction  of  h, 
computed  at  the  point  Xq.  The  Gateaux  derivative  is  not  necessarily 
continuous  with  respect  to  h either  in  the  strong,  or  even  in  the 
weak  topology  of  H2"  It  only  implies  the  existence  of  a linear 
operator  ^ ^ such  that  for  a fixed  h ^ H2  and  for  sufficiently 
small  (in  absolute  value)  t g R, 


f(xQ  + th)  - f(xQ)  =tL^^  jh  + ‘Z(xQ.th) 


(1.8) 


where 


lim 

t-^-0 


»j(Xo;th) 


0. 


(or  weak  lim. 
t-*0 


-£  't ( Xq ; th ) = 0. ) 


We  would  like  to  point  out  that 


lim 

t-»0 


7.(  Xpfth) 

t 


0 


lim 


IlhIkO 


||7(  XQ;h)|| 

( m 


does  not  imply 


0 


for  arbitrary  h € H2  (of  suf f icieiitly  small  norm)  . In  other  words 

the  Gateaux  derivative  may  exist  in  the  direction  of  a vector  hj^ , 

but  may  fail  to  exist  in  the  direction  of  a vector  h2*  This  is 

easily  checked  to  be  true  even  in  some  two  dimensional  cases) 

2 2 

f : R 


1 . 4 Frechet  differentiation • 

Assume  that  the  operator  L in  (1.8)  is  linear  and  bounded,  i.e. 
for  fixed  Xq  € H2  we  have 

i)  L (oh,  + gh  ) = ctL  h^  + gL  h„ , and  there  exists  M > 0 


X 1 P X 7 ' 
Xq  J.  ^0  ^ 


such  that 


ii)  1|l  hll  < M(x«)  !lh|!  for  all  h of  sufficiently  small  morm. 

Xq  - U 

We  also  assume  that  for  all  such  h g ^2  true  that 

iii)  f(xQ+h)  - f(xQ>  = h + where 


IlhTf-vO 


"«(xQ;h) 


= 0, 


( Jee  C 7 J J 


(1.8) 


Then  L is  called  the  Frechet  derivative  of  the  operator  f eval- 
^0  . 
uated  at  Xq.  The  following  condition  can  replace  condition  (iii) 

^(xQ,th) 

iv)  lim  = 0 uniformly  with  respect  to  h on  all 

t-»0  t 

bounded  subsets  of  H2  • 

The  fact  that  (iii)  and  (iv)  turn  out  to  be  equivalent  is  not  trivial. 


The  proof  of  it  can  be  found  in  [21 ]. 


Theorem.  Frechet  derivative  exists  if  and  only  if  the  Gateaux 
derivative  exists  and  satisfies  condition  (iii).  The  proof  is  easy 


8 


and  will  not  be  given  here.  (See  I52 ] • ) 


1 . 5 Fr^chet  differentiation  of  functionals . 

In  this  case  we  consider  the  Frechet  derivative  of  the  map 
f:  c H.  -+•  R.  Let  the  Frd’chet  derivative  L exist  at 

1 Xq  u 

Then  L is  a continuous  linear  map  from  fi  ^ H.  into  p.  Therefore 

for  some  h € H.  we  have  L ch  = c L h,  V c € R,  and  by  the  Hahn 

Banach  theorem  L can  be  extended  to  all  of  H,  without  changing 

Xq  1 

its  norm. 

By  the  Riesz  representation  theorem  there  exists  a unique 
Z € H.  such  that  L h = <Z,h>  for  all  h € h, . Hence  for  all  h € H 

± Xq  X 

the  Fre'chet  derivative  of  a functional  f evaluated  at  Xq  is  a 

continuous  linear  functional  L , and  L h = <Z  ,h>  (1.9) 

^0  ^0  ^0 


for  all  h ^ Hj^ . 

The  unique  vector  Z is  called  the  gradient  of  f evaluated 
at  Xq.  Clearly  Z depends  on  the  definition  of  the  inner  product 
in  Hj^ . However  the  existence  of  the  gradient  depends  only  on  the 
existence  of  the  Frechet  derivative  which  in  turn  depends  on  con- 
vergence of  our  limiting  process.  Convergence  is  a topological 
property,  and  remains  invariant  in  equivalent  topologies.  Hence 
the  introduction  of  an  equivalent  norm  in  Hj^  will  not  affect  the 
existence,  or  non-existence  of  a gradient  (or  of  the  Frechet  deriva- 
tive). In  fact  let  os  introduce  a new  inner  product  (the  so  called 
energy  product,  and  corresponding  energy  norm)  [ , ] 

[x,y]  = <Tx,y>,  where  T is  positive  definite, symmetric  operator. 
Denoting  by  |||  xj|l  = [x,x]  = <Tx,x>,  and  by  ||y||  = <y,y>,  we 

A 4 

have  <grad  f,  h>  = [grad’^f,h]  = <Tgrad*f,h> 


(1.10) 
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» 

where  grad'  is  the  gradient  of  f in  the  topology  induced  by  the 

norm  |j|  ||| . Since  relation  (1.10)  is  true  for  all  h € we  have 

# -1  # -1 
the  equality  T*grad  = grad.  If  T exists,  then  grad  = T *grad. 

If  the  operator  T is  bounded  away  from  zero  (i.e.  the  spectrum  of 
T does  not  have  zero  as  a limit  point)  then  there  exists  a constant 
Y > 0,  such  that  for  any  x ^ n .^(T)  ||  x||  ||j  x ||| . (-^  (T)  denotes 

the  domain  of  T) . Hence  convergence  (or  continuity,  or  existence 
of  limit  points)  in  the  new  norm  1|1  j!|  implies  convergence  (or  con- 
tinuity...) in  the  old  norm  ||  ]]  . 

In  the  final  dimensional  case  all  these  statements  are  trivial, 
and  all  norms  <,>,  [ ] where  <Tx,y>  = [x,y] , and  T is  positive 


definite,  and  symmetric,  are  equivalent  and  generate  the  same 
topology. 

At  this  point  it  is  appropriate  to  make  a comment  about  the 
non-existence  of  the  Frdbhet  derivative  in  some  Hilbert  space 
Suppose  the  Gateaux  derivative  of  $ (in  arbitrary  direction  h)  exists, 
but  fails  to  be  a continuous  functional  of  h in  the  usual  (norm) 
topology  of  i.e.  we  can  find  a sequence  of  vectors  {h^ } € 
such  that  lim  H h - h^^!!  = 0,  but  lim  $(hj^)  ^ $(h)  . It  may  be 

i.>-oo  i^voo 

possible  to  introduce  a different,  and  non-equivalent  inner  product 

{,},  with  a corresponding  norm  !!  11(2)  that  |i(h)  is  a continuous 

functional  of  h in  the  new  (norm)  topology. 

If  the  inner  product  is  fixed,  and  no  attempt  will  be  made  to 

change  the  topology,  then  we  can  use  the  sloppy  notation 
54-(x) 

or  — r — ^ to  denote  the  gradient  of  ! evaluated  at  x«,  where 

6 X Xq  V 

$ is  a Frechet  differentiable  functional  whose  values  depend  on 
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X e H,  , i.e.  H,  R.  L (x)  = <Z  ,x>  = I ,x> 

X X Xq  Xq  6X  'Xq 

is  the  corresponding  Frechet  functional.  This  notation  is  both 
confusing,  and  sloppy,  but  it  has  been  consistently  in  use,  and  the 
author  of  this  paper  is  also  guilty  of  having  abused  it. 

It  is  partially  vindicated  by  a physical  explanation.  The 

Hilbert  space  setting  has  obscured  our  vision  to  a certain  extent, 

since  we  identify  H with  its  dual,  and  we  are  allowed  to  neglect 

the  bookkeeping  which  vector  is  in  H,  and  which  vector  is  in  H* 

(the  dual  of  H,  which  is  identified  with  H) . This  distinction 

becomes  inportant  later  on  when  we  try  to  interpret  certain  equations 

of  mathematical  physics.  In  the  remainder  of  this  paper  we  shall 

use  the  notation  grad  $ if  $-is  a functional  $:  ->  R.  We  shall 

also  find  examples  where  $ is  a functional  $:  ® K2  ®.  . .©  R 

(or  <!),.  that  is  $ maps  the  direct  product  of  Hilbert  spaces  H^^ , K2  • . . 

Hj^  into  the  real  (or  complex)  numbers . In  this  case  the  notation 

grad  . is  too  cumbersome  and  we  shall  write  J to  denote 

vn.; 

f(x-+th)  - f(x-) 

L (h)  = = <1^"*  ( f 

X X t+0  ^ 

with  h € H^,  Xq  € H^^  ® H2  ®...®  Hj^.  That  is  we  generalize  simultan- 
eously the  idea  of  a gradient  and  of  a partial  derivative,  by  allow- 
ing the  variation  of  f to  take  place  exclusively  in  H^,  and  regard- 
ing the  components  of  x in  Hj , j ^ i,  as  fixed.  We  shall  refer  to 
this  derivative  as  gradient  of  restricted  to  the  space  H^ . This 
idea  plays  a cjOE^rial  part  in  formulation  of  multiple  variational 
principles.  For  exanple  the  following  equations  of  mildly  nonlinear 
elasticity  corresponds  to  critical  points  of  a simple  functional 


11 


f 


restricted  to  respective  subspaces  of  the  space  H = 7 H. . The 

e ^ 

Maxwell  and  Morrera's  stress  function  relations,  the  equations  of 
equilibrium,  the  equations  of  compatibility,  Ricci's  equations. 
(See  [o'3]).A11  become  easily  interpreted  in  this  context. 


1 . 6 Some  rules  of  manipulation  for  gradients  of  functionals 
We  shall  restrict  ourselves  first  to  the  simplest  case:  41:  H R, 

^ H,  4>  is  Frechet  differentiable  in  ^ H. 

We  shall  denote  by  4-^^^  the  gradient  of  i = Z € where 
<Z,n>,  Z € the  Frechet  derivative  of  <5  coinciding  with 

the  Gateaux  derivative  in  the  direction  n»  We  conpute  the  value 
of  the  gradient  at  a point  x^  €.5.  = Z(Xq).  If  Z(Xq)  = 0 

(the  zero  vector)  we  say  that  Xq  is  a critical  point  of  the  func- 
tional 4>.  A basic  theorem  of  optimization  theory  (Vainber^  (49^) 
states  a necessary  condition  for  stationary  behavior  of  a functional 
<I>  (particularly  for  a maximum  or  minimum)  at  the  point  Xq  € H.  The 
point  Xq  must  be  a critical  point  of  J,  if  5 is  Frechet  differenti- 
able at  Xq.  In  particular  this  is  a necessary  condition  for  a 
local  maximum  or  minimum  of  a Frechet  differentiable  functional. 


Example  1 . $ = <Yj^ ' ^2^ 

^2  € H,  X , 

^2  " ^2  Z = 4^  = (C^  + C^)  X. 

Proof  of  this  statement  follows  from  the  definition.  We  compute 
the  Gateaux  derivative  of  ? in  the  direction  of  a vector  n ^ H. 


= (C^  + C2)  <n,x>  = <n/  (Cj^  + C2>X>. 


Hence  4>jj  = + C2)X. 

Example  2.  4i(X)  = * *2^^^  Then  the  gradient  of  4>: 


(1.11) 
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This  formula  can  be  obviously  generalized  to  polygonal  paths.  The 
common  way  to  introduce  integrals  of  vector  fields  along  an  arc 

3 

re  R is  to  consider  a subdivision  of  the  arc  Pr  1. 1 by  points 

ia,bj 

Pq  r P2‘**'Pn  introducing  the  Riemann  sums  ; 


n 


.1 

i=0 


f(  q)  , 


where  is  a point  contained  on  the  interval  [p^^,  ^i+1^* 

3 3 

f;  R ->  R is  a continuous  map  and  if  r is<  "reasonably  smooth". 

(For  example  we  can  allow  a finite  number  of  corner  points  without 
causing  any  problems) . Such  Riemann  sums  generalize  to  Hilbert 
(Banach)  spaces  by  simply  replacing  the  three-dimensional  inner 
product,  by  the  appropriate  inner  product.  If  f is  a continuous 
map  f:  .H  H,  and  r is  a smooth  arc  we  can  easily  show  the  conver- 
gence of  the  Riemann  sums : 


lim  { T <(Pi+i  - Pi»  f(Ci)>).  (1.12) 

IIPi+2^-Piii--0  i=0 


i In  the  case  when  r is  a parametrized  arc  we  can  choose  the  points 

V to  coincide  with  p.,  and  the  integral  (1.12)  can  be  reduced  to 

r the  form 

1 

/ <t)(t)dt. 

0 

i' 

1.8  A basic  condition  for  the  existence  of  a potential. 
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i 


i 


% 


t 


a region  of  H,  such  that  for  any  € n the  Integral 

(Pl»P2) 

j <f(x),dx>  = j (|i(t)dt 
r - ‘ ‘ t{p^J=t^ 

is  independent  of  r,  for  any  parametrized  arc  r,  and  only  depends 
on  the  end  points  pj^,  p2«  Then  there  exists  a functional  $(  x)  such 
that  f(x)  is  the  gradient  of  J. 

Proof . Choose  a (fixed)  point  Xq  € ^ H.  Define  for  p € jj  the 

functional 

1 

4>(p)  = /.<f(xQ  + t(p  - Xq)),  (p-XQ)>dt  + 4>(Xq), 

and  (without  any  loss  of  generality)  put  '^(Xq)  = 0. 

If  fi  contains  the  origin,  we  may  choose  Xq  = 0 , and  define 
1 

♦ (x)  = / <f(tx),  x>dt. 

0 ~ ~ 

It  is  easy  to  check  that  f(x)  is  the  gradient  of 

Note:  Observe  that  we  did  not  make  the  assumption  that  f:  H -*•  H 

is  linear. 

There  is  the  obvious  problem  of  knowing  how  to  apply  this  con- 
dition. How  do  you  check  that  this  integral  is  independent  of  the 
path  for  every  path  in  r?  In  most  cases  this  turns  out  to  be  an 
impossible  task.  At  this  point  we  recall  another  condition  for  the 

existence  of  a potential  function  in  Assume  that  Q is  simply 

3 3 

connected,  that  f:  R R is  smooth  and  curl  ( f ) exists  at  all 

points  of  fi,  and  curl(fjs  0 in  n.  Then  f is  a gradient.  Of 
course  we  do  not  know  what  curl (f) means  in  an  infinite  dimensional 
space.  However  in  R curl(f)H  0 if 


1 
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f 1 — Ir2r3i 


This  condition  generalizes  quite  easily  to  infinite  dimensional 
space.  However  we  shall  have  to  delay  its  discussion,  since  we  have 
not  defined  yet  the  meaning  of  second  (or  higher)  order  derivatives 
of  a functional.  Observe  that: 


3 3 X j 


3Xj3x^  * 


3 2 

This  condition  is  trivially  satisfied  in  R if  4>  € C (fi),  (this  is 
known  as  Tonelli's  theorem).  In  a Hilbert  space  there  is  no  obvious 
reason  why  second  derivatives  should  commute,  and  this  condition 
assumes  a much  deeper  meaning. 

The  generalizations  of  Cauchy's  integral  condition  in  higher 
dimensional  complex  spaces  are  too  involved  to  be  of  practical  use 
in  Hilbert  spaces.  At  least  this  is  the  opinion  of  the  author, 
which  may  of  course  turn  out  to  be  wrong. 


mccsDnc  pj 


Chapter  2 


2.1  The  Euler-Lagrange  equations , and  critical  points  of  a functional. 


The  Euler-Lagrange  equations  for  the  problem  of  extremizing  the 


value  of 


/ ^ x,x' , t)dt 


8 ^ _ d I 8^  \ ^ 

8X  dt  ' 3x'  ' “ 


(2.1) 


The  Legendre  transformation  introduces  tha  generalized  monentum 


p = — . 
^ 3x' 


(2.2) 


Hence  (2.1)  can  be  written  as 


8id  _ • _ _ 8H 
8X  ^ ” 8X 


(2.3^) 


= X 


where  H = x p - 


(2.3^) 


We  are  now  ready  to  formulate  the  Hilbert  space  analogs  of  equations 
(2.3^),  (2.3^),  and  (2.1).  The  basic  ideas  of  generalization  of 
equations  (2.3^),  (2,3^)  go  back  to  Korn  who  observed  dual  varia- 
tional principles  in  theory  of  elasticity  of  materials  not  necessar- 
ily obeying  Hooke's  law.  Friedricks  has  interpreted  Korn's  results 
in  terms  of  modern  operator  theory  13  1. 
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An  intuitive  approach  to  the  postulation  of  duality  and  formu- 
lation of  dual  variational  principles  originated  with  Noble  [35 ] , 
and  has  been  successfully  applied  by  a number  of  authors  to  problems 
of  mathematical  physics  (see  for  exanple  the  monograph  of  A.  M. 
Arthurs  13]). 

2 . 2 Generalized  Solutions,  energy  norms  and  extremal  points  of 
functionals . 

rl,  ' * ' 

The  well  known  case  of  a close  relation  between  a generalized  solu- 
tion o equation: 

Au  - f = J2f  (2.4) 

in  a Hilbert  space  and  extreme  value  of  a functional^ when 

A is  a positive  definite  operator,  bounded  away  from  zero,  and  the 

domain  of  A is  dense  in  H.  In  this  case  we  can  introduce  a new 

2 

inner  product  [u,v]  = <Au,v>,  and  a new  norm  !11 111  ~ <Au,u>  . 

Let  denote  the  closure  of  the  new  inner  product  space.  We 
excimine  the  corresponding  functional: 

P(u)  = <Au,u>  - 2<u,f>  = ||j  u “ 2<u,f>  (2.5j 

2 2 

Since  A is  bounded  away  from  zero  (meaning  <Au,u>  |1  u||  for 

ff  fll 

some  C > 0 , or  l||u|||  > C ||u]|)  it  follows  that  |<u,f>|  5. ill  u Hj  . 
Hence  by  Riesz  representation  theorem  there  exists  Uq  € H^  such 
that  <u,f>  = [u,Uq],  and  F(u)  = |||  u - Uq  ||1  ^ - jll  Uq  11]^ 

Hence  F(u)  attains  a minimum  at  Uq,  and 

min  F(u)  = - |||  Uq  |H  ^ 
u CH^ 


(2.6) 
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As  before  denotes  the  Hilbert  space  obtained  by  closure  of  the 
inner  product  space  with  the  product  [•,*].  The  point  Uq  € is 
called  the  generalized  solution  of  the  equation  (2.4).  If  Uq  was 
in  the  domain  of  A,  then  it  is  a "genuine  "(or  classical)  solution. 
Otherwise  it  is  called  a generalized  solution  of  (2.4).  In  fact 
since  A is  positive  definite,  there  exists  an  operator  B,  such 
that  A = B*B,  and  for  all  u € u € and  <Au,u>  = <Bu,Bu>. 

We  have  the  containirent  ^ • However  ^ As  does  not  wean 

A c U o 

that  Uq  is  necessarily  in  the  domain  of  A.  As  an  exaitple  consider 


A = - , .»  = C^[0,1]  C L^[0,1],  u(0)  = u(l)  = 0, 

dx 


The  generalized  solution  of  5 u = f corresponds  to  the  mini- 

dx'^ 

mization  problem  for  the  (L2)  functional 


<-  — 5 u,  u>  - 2 <f,u> 

dx 


^ (u'*  • u)dx  - 2 ^ (fu)dx 


^ (u')^dx  - 2 (fu)dx. 


The  problem  is  that  the  existence  of  the  functional 


F(u)  = 


•)^dx  - 2 / (fu)dx 


requires  only  u'  € L2l0,l]  and  does  not  require  twice  differentiability 


- 
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B = i is  one  of  the  possible  corresponding  square  roots  of  the 

d ^ 

operation  of  the  operator  - — 5^.  Hence  the  minimization  problem 

dx'^ 

can  be  accomplished  in  the  Sobolev  space  H^[0,1]  instead  of 

u € C [0,1]  c:  L2(0,1].  The  energy  space  turns  out  in  this  case  j 

to  be  H^[0,1].  Of  course  in  the  above  case  we  had  a lot  of  pro- 
perties of  the  operator  A,  (positive  definite,  bounded  away  from 
/ 

zero)  which  made  life  very  easy,  and. gave  us  so  easily  an  equiva- 
lent extremal  problem  for  the  functional  F(u) . A feature  of  this 
functional  which  we  have  not  mentioned  yet  was  convexity. 

Definition.  Let  <t:  n = H -v  R be  defined  on  a convex  set  n = H. 

Then  <Ji  is  called  convex  if  for  any  Xj^ , X2  € the  following  in- 
equality is  true: 


V 

J. 


4>(Ax^  + (1-A)x2)  A't(x^)  + (1-A)$(x2)  0 < A < 1. 


If  strict  inequality  holds,  then  f is  called  strictly  convex. 
li  is  called  concave  if  -'J  is  convex.  We  observe  that  any  norm  is 
a convex  functional,  since  it  obeys  the  triangular  inequality. 


► 

r*: 


t 


Hax  + (l-A)yl!  <_  a|!x||  + (l-A)llyll. 

A study  of  convex  functionals  can  be  reduced  to  a study  of  convex 
sets  in  H,  by  introduction  of  the  "epigraph"  of  the  fujc^tional  $, 
that  is  of  a subset  of  R x H 

[4>,C]  = {(r,x)  €r  X H:  4-(x)  < r,  x *-H} 
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where  C is  a convex  subset  of  H.  (see  Luenberger  [55]  for  an 
exposition) . 

2 . 3 Conjugate  convex  functionals . 

Let  be  a convex  functional  4>:  n s B ->  R,  defined  on  a convex 
set  n c B-,  where  B is  a Banach  space  and  B*  is  its  dual.  The  set 
* conjugate  to  n relative  to  $ , B*,  is  defined  as 

n*  = {p  6 B*:  sup  <p,x;  - $(  x)  < m for  some  M € R} 

' x€n 

The  conjugate  functional  4i*(p)  is  defined  by  the  relation 

<l>*(p)  = sup  (<x,p>  - 4>(  x)  } . 

x€0 

Proposition.  If  is  convex,  then  so  is  $* . Proof . For  any 
number  A,  0 < A < 1,  check  that 

sup  Ap,  + (1-A)P5>  - $(  x)  } 

x€n 

£ A sup  {<x,  p,  > - <!>(  x)  } + (1-a)  sup  {<x,  p_>  - $(x)}  . 
x€n  x6a 

A geometric  interpretation  of  the  convex  functional  it*  can 
be  intuitively  argued  as  follows. Consider  a convex  set  n c p x b, 
where  B is  a Banach  space.  A family  of  hyperplanes  in  R « b is 
given  by  an  element  of  R x B*,  i.e.  an  ordered  pair  (s,  p)  S 6 R, 
p €b*.  a particular  hyperplane  is  obtained  by  setting 
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<(r,x),  (s,p)>  = c,  r € R,  x € B,  or  rs  + <x,p>  = c,  where  c is 
some  real  number.  By  rescaling  we  can  always  choose  s = -1,  pro- 
vided s 0 , Then  C = <x,p>  - r defines  a family  of  hyperplanes 
of  R X B.  4>*(p)  = sup  (C)  is  the  sup.  of  values  for  which  't*(p)  is 
a support  hyperplane  of  [ '!>,  C]  . 

Convexity,  and  concavity  of  certain  functionals  are  closely 
tied  with  existence  theory  of  maxima,  and  minima.  For  the  time 
being  we  shall  only  concentrate  on  sufficiency  conditions  which 
generalize  the  Euler-Lagrange  equations  of  classical  calculus  of 
variations . 

2 . 4 Critical  points  of  fucntionals  and  the  equation  Ax  = f . 

We  shall  consider  a fairly  arbitrary  linear  operator  A mapping 
.5^  c into  where  is  a Hilbert  space.  We  wish  to  solve 

(for  X € H^)  the  equation  Ax  = f (f  C , where  we  shall  assume 
that  A = T*T  (A  is  positive  definite)  and  the  domain  of  T is  dense 
in  Therefore  T*  (the  adjoint  of  T)  is  uniquely  defined.  T, 

T*  are  the  linear  maps : 

T;  H.  H2.  T*;  H2  -»• 

Hence  the  equation  Ax  = f (€  H^^)  can  be  rewritten  as  a pair  of  equa- 
tions 

T X p (€  H2)  , (2.7) 


T*  p 


f 


(c  Hj^). 


(2.8) 
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We  consider  the  Hilbert  space  H = ® H2  whose  elements  h are 

ordered  pairs  h = (^»P)  p € ^2)  t with  the  product 

{hi,  ^2  } = (x^,  x^)  + <pj^,  P2> 

where  (,)  is  the  inner  product  in  , and  <,  > is  the  inner  product 
in  H2.  Consider  arbitrary  vectors  p € H2 , x € i.e.  an  arbitrary 
w = { x,p  }6  H and  the  corresponding  value  of  the  functional 
L:  H -)■  R 


L(x,p)  = <Tx,p>  - <p,p>  - (f,x) 


(2.9) 


It  is  a fairly  trivial  result  that  provided  the  gradient  is 
uniquely ' defined  at  the  point  Wq  = { Xq , Pq}  € H,  then  the  equations 
(2.7)  and  (2.8)  are  necessary  and  sufficient  conditions  for  vanish- 
ing of  Hence  the  functional  L has  a critical  point  at  Wq  if 

and  only  if  the  equations  (2.7)  and  (2.8)  are  satisfied. 

Again  the  condition  x € has  been  replaced  by  condition 
X € which  is  in  general  easier  to  satisfy  since  A = T*T.  (See 

the  next  section  for  a more  rigorous  statement.)  A more  detailed 
look  at  equations  (2.7)  and  (2.8)  show  that  if  we  fix  p € H2  and 
vary  only  x € H^,  = 0 (in  H^)  if  and  only  if  (2.8)  is  satisfied, 

and  vice  versa  if  we  fix  x € H^^  and  vary  p € H2 , (in  H2) 

if  and  only  if  (2.7)  is  satisfied. 

These  observations  have  been  originally  implied  by  Kato  [22] • 

(We  are  using  the  same  notation  denoting  the  gradient  of  the 


functional  L(x)  (with  p a fixed  vector  of  ^2^'  always  assuming  that 
this  gradient  exists,  (in  a fixed  topology  of  the  space  I ) and  if 

= 0 at  X = Xq  then  is  defined  in  some  neighborhood  of  the 

point  Xq  in  the  appropriate  topology.  Unless  otherwise  specified  we 
shall  use  the  topology  induced  by  the  norm.  Again  we  observe  that 
Wq  = {Xq,Pq}  which  is  a critical  point  of  the  functional  L does  not 
have  to  be  a "genuine  solution  of  the  original  equation  Ax  = f, 
since  x»  may  be  in  the  complement  of  the  domain  of  A. 

The  system  (2.7),  (2.8)  has  been  designated  by  Noble  ([35]*^)  and 
Rail  ( ] ) as  a Hamiltonian  system.  The  name  is  easily  explained 
if  we  call  the  functional  j + <f,x>  = W(p,x)  the  Hamiltonian 

and  observe  that  (2.7),  (2.8)  can  be  rewritten  as 


T X = W 

P 

(2.10) 

T*p  = 

(2.11) 

which  is  the  Hamilton's  system  of  canonical  equations  in  the  special 
d d 

case  ~ ~ ^ * with  x being  the  vector  of  generalized 

displacements,  and  p of  generalized  momenta. 

In  anology  with  the  terminology  of  classical  mf:chanics  the 
fucntional  L will  be  called  the  action  functional  cx  the  Lagrangian 
(integral)  functional.  We  observe  the  following  peculiarity  of 
our  discussion:  The  equation  Ax  = f is  an  equation  in  . The 

introduction  of  the  "splitting  space"  II2  was  a result  of  the  decom- 
position of  the  operator  A;  into  T*T:  T:  ^ K2,  T* : H2 
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In  particular  the  work  of  Browder  and  Gupta  indicates  that  under 
certain  conditions  while  is  a Banach  space,  H2  may  be  chosen 
to  be  a Hilbert  space,  with  resulting  theoretical  advantages.  In 
this  paper  will  be  assumed  to  be  a Hilbert  space.  To  make 
rigorous  the  ideas  loosely  expressed  here  we  need  some  additional 
definitions  from  functional  analysis. 

Definition  Let  T be  a (linear)  operator  mapping  some  subset  of 
a Banach  space  Bj,  into  a Banach  space  B2 . The  graph  of  T,  denoted 
by  G(T)  is  a subset  of  B^^  x B2  consisting  of  all  ordered  pairs  of 
the  form  {u  ,Tu),  u € .5^,  The  following  statement  is  known  as  the 
closed  graph  theorem.  Let  = B^.  Then  T is  bounded  if  and  only 
if  the  graph  of  T is  closed. 

Definition  (Kato  [5^])  Any  linear  mainfold  C contained  in  ^ B^^ 
is  called  the  core  of  T if  the  set  {u,  Tu},  u € c,  is  a dense  sub- 
set of  G(T).  An  important  theorem  due  to  Von  Neumann  asserts  the 
following.  Let  T be  a closed  operator  T:  ->•  H2;  dense  in 

Then  A = T*T;  -+  is  a self  ad  joint  operator  and  is  a 

core  of  T.  Hence  D^  2 if  the  containment  is  proper  i.e. 

D.J,  ^ D^  then  the  corresponding  variational  solutions  describing  the 
original  problem  in  terms  of  critical  points  of  some  functionals 
in  H2  will  exhibit  solutions  which  are  not  "genuine"  solutions, 
that  is  they  are  not  in  the  domain  of  A. 

2 . 5 The  Legendre  transformation 

In  classical  mechanics  we  encounter  the  following  transformation: 


' and  generalizing  this  concept  to  an  abstract  Hilbert  space  setting. 

1 

Again  perhaps  the  best. starting  point  is  the  equation 

I Ax  = f ( € Hj^)  , A: 

i 

j ' 

for  which  we  intend  to  establish  (weak)  solutions  which  correspond 
t j to  a critical  point  of  some  functional  L,  or  to  two  (or  more) 

critical  points  of  functionals  in  some  spaces  possibly  other  than 


The  case  we  have  discussed  already  presumes  that  A is  positive 
definite  and  that  we  can  find  a soace  H..  such  that  x € Tx  = o 
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Of  course  the  choice  of  is  non-unique,  as  can  be  seen  by  studying 
even  the  simpl  est  examples 

2 2 
(-  ^ : n^(R)  ->  L,(R),  - ^ = T*T, 

dx^  dx^ 

T = i ^ , T*  = i ^ , T;  H^(R)  H^(C)  , 

T*:  Il^I)  - Lj(R)  , or  T = |j  , T*  = - . 

T:  H^(R)  H^(R)  , T*:  H^(R)  ->■  , 

showing  that  A Can  be  factored  thru  different  Hilbert  spaces 

and  . More  involved  exeunples  will  be  offered  later  of  multiple 
factoring  process  with  corresponding  multiple  variational  principles, 
in  which  the  choices  of  intermediate  Hilbert  spaces  are  not  uniquely 
determined. 

We  shall  now  consider  the  effect  of  the  vector  Tx  € H^^  on  the 
value  of  the  Lagrangian  functional  L.  Regarding  p as  a fixed 
vector  in  H^^,  we  consider  L(x,y)  as  functional,  mapping  the  pair 
X € y ^ H2  into  the  real  line,  where  y = Tx. 

Ly  denotes  the  gradient  of  L (whenever  it  exists)  restricted 
to  H2*  (x  € is  ignored).  Similarly  denotes  the  gradient  of 
L in  It  is  a straightforward  computation  that  with  p regarded 

as  fixed,  and  with  y = Tx,  that 


V 


(2.12) 
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and  that 


<Tx,  p>  - L = ^ <P»P>  + (f»x)  ^ W (2.13) 

defines  a new  functional  W(x,p)  satisfying 

N 

(2-/4) 

The  notation  is  the  same  as  before.  is  the  gradient  of  W restric- 
ted to  the  space  and  is  the  gradient  of  restricted  to  the 

space  H2 . We  s)iall  describe  the  relation?  (2.12),  (2.13)  as  the 
Legendre  transformation.  /'Jee  <'  ciosSi<:.e.i  ^ 


T 


w 


Chapter  3 

3.1  An  example  of  multiple  critical  points,  and  applications  to 
the  general  theory  of  solids. 

We  shall  consider  the  equations  of  classical  elasticity,  assuming 
small  strain  theory,  but  not  necessarily  small  rotations.  The 
equations  of  equilibrium  assume  the  form: 

ir  "xx  ^®xy  - "xy  ^®xz  ^xz^ 

+ ®xx^  ^xy  ^®xy  ” “z^  Vy  ^ 

+ 1?  ®xx^  "xz  <®xy  ■ “z^  ^yz  ^®xz  “y^  ^zz^  = ° 

(3. 


Equations  (3.2),  (3.3)  are  obtained  by  cyclic  permutation  of  letters 

X,  y,  z.  Here  is  the  linear  strain  matrix  e^x  ^ ®xx  “ ^ ' 

^ :e  = u.  are  the  rotation  components 

^xy  xy  2 6y  3x  ' i 


i / aw  _ iZ  \ 

'“x  “ 2 ' 3y  3z 


1 , 3U  _ 3V  t 
’y  “ 2 ' 3Z  3X  ' 


_ 1 , 9V  _ 3U  , 
“z  2 ' 3X  3y  ^ 


(3.2) 


i .e. 


-►  1 
w = -o 
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If  we  denote  by  o the  following  stress  tensor 
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we  can  formulate  the  following  sets  of  equations  of  equilibrium 


3x 

0 

0 

3 

ay 

0 

3 ■ 
3Z 

' “'xx  ' 

B a = 

0 

3 

ay 

0 

3 

3X 

3 

3Z 

0 

“yy 

0 

0 

3_ 

3Z 

0 

3 

ay 

3 

3X 

“zz 

= 0, 


(3 


1 

32 

However  the  pairs  {U,4>},  {o,e>,  {K,X}  are  not  vectors  (in  the  app- 
ropriate spaces)  which  can  be  picked  independently  of  each  other. 

If  U denotes  the  displacement  vector  B*U  = e is  the  definition  of 
the  linear  strain. 

B A*B*U  = 0 (3.12) 

is  exactly  the  set  of  Ricci's  equations.  (see  Washizu  [20] 
for  explanation  of  their  importance  in  linear  elasticity.)  The 
following  diagram  illustrates  the  relations  (3.5)  - (3.12) 


e 


iil 


f\J  € H3 


I 

If 


They  are  related  to  each  other  through  the  constitutive  equations 
of  the  solid.  For  excunple  if  Hooke's  law  is  assumed,  then  we  have 


9 X 9 matrix  of 


anisotropic  coefficients) , 


‘'ij  ■ ^ijkH  “k4  ' 


or 
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where  is  a constant  9x9  matrix  such  that  ^ 

(the  9x9  identity  matrix) 

A well  known  set  of  variational  principles  in  classical 
elasticity  can  be  derived  by 

a)  Assumming  a constitutive  equation  of  a solid 

b)  Restating  all  basic  equations  of  elasticity  as  the  existence 
of  a critical  point  of  the  corresponding  functionals: 


= W - 


= W - 


2 ^ 


2 ^ 


♦,  = w - 


♦,  = w - 


2 < Ax,e  > 


2 ^ ^ (H,) 


■-) 


(3.13) 


1. 


= W - < U,Ba  > 


(H3) 


= W - < B*U,a  > 


(H2) 


♦-  = w - 


< U,q>  > 


■ 

P. 

i'm 


A. 


Of  course  we  have  deliberately  ignored  the  basic  problem  of  solving 
the  basic  boundary  value  problems  of  Elasticity  (using 
Muskhelishvili ' s terminology)  and  concentrated  on  the  formulation 
of  fundeunental  principles  with  natural  boundary  conditions. 

(See  definition  and  discussion  in  section  a / ) 

In  fact  should  another  loop  be  added  to  our  mapping  diagram 


K 

t. 


(3.1)  we  can  immediately  formulate  the  corresponding  critical  point 
statement,  hence  a variational  principle  for  the  corresponding 
functional. 


3.2  Possible  (future)  applications  to  the  basic  theories  of 
solid  materials. 

The  basic  thermodyncunics  laws  applicable  to  solid  materials  can 
be  summarized  in  the  equations  (3.13),  (3.14)  and  inequalities 
(3.15),  (3.16)  below.  The  solid  occupies  a region  E^,  in 
which  the  following  relations  are  valid. 

(3.14)  E = I (Vp)T^j  + Q , 

(first  law  of  thermodynamics) 


(3.14®) 

(3.14^) 


4>  = iMC^^,0,  q“)  , 


T.  . = T..  (c„„»0»  q“) 

1]  1]  mn  ^ 


(the  principle  of  material  indifference) 


0 > 0 , 
0S  > 0 


(Absolute  temperature  is  non-negative) . 
(Positive  rate  of  dissipation) . 


(3.15) 


h^0,^  0 (3.16) 

(Heat  conduction  in  the  direction  of  negative  temperature  gradient) . 
The  symbols  used  here  have  the  following  meaning.  E is  the  internal 
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energy,  C^j  is  the  Cauchy-Green  strain  tensor,  p is  the  mass  density 
p(0)  = pQ  / ® temperature,  S entropy,  \p  the  free  energy,  h^  the 
heat  flow,  Q the  heat  supply  per  unit  of  mass. 

Dots  denote  differentiation  with  respect  to  time,  commas  - 
covariant  derivatives.  Summation  notation  is  used  unless  otherwise 
stated,  are  additional  independent  variables,  called  internal 

variables,  which  are  not  necessarily  observable,  i.e.  the  implicit 
variable  theorem  may  be  ''^applicable  to  relations  (3.14). 

It  has  been  shown  by  Velanis  I4i 1 and  others  that  the  stress 
distribution  consistent  with  the  first  and  second  law  of 
themodynamics  satisfies  the  equation 
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i.e.  T . . is  the  gradient  of  ^ W in  the  appropriate  product  of 
Sobolev  spaces,  which  are  subsets  of  L2(fi)  L2  (jT)  . Recalling 
that  the  entropy  S satisfies  the  relation 


- S 


M 
86  ' 


(3.17) 


and  that  the  temparture  0 is  always  positive,  the  first  law  of 
thermodyncunics  can  be  written  as : 


-S 


dt 


86 


- lilL  4“ 

8q“ 


(3.18) 


(which  makes  no  reference  to  constitutive  equations  for  the  material.) 
We  shall  attempt  to  follow  at  this  point  a fascinating  idea  of 
Ilyushin  [1^],  who  suggested  that  the  material  has  no  concept  of 
our  idea  of  time,  and  that  processes  within  a solid  should  be 
parametrized  with  respect  to  a "material  time",  which  was  later 
defined  by  Riviin  by  the  relation  dr  = (dC^^  (I40l). 

Also  see  [4y]  , [421*  particular  the  internal  variables  q 
should  only  depend  on  t,  and  preferably  in  a simple  manner. 

We  assume 


T = q“(t)  , •|(t)  ,t) 


(3.19) 


8t 


HI 


8C 


iD 


8q 


8t 

a 


8t  r 
8^ 


8t 

8t 


and 


dx 

at 


(}-70) 


We  hypothesize  that 


> 0 , (3.21) 

( 

which  is  a sensible  physical  assumption,  since  otherwise  the 

material  would  react  to  future  physical  condition,  or  could,  in 

a manner  of  speaking  predict  the  future.  The  other  assvunption 

frequently  made  here  is  the  existence  of  a steady  state  condition. 

The  assumption  that  C = 0 is  possible  independently  of  q®,S,6, 

or  t,  or  q“  = 0,  etc.  implies  that  each  term  could  be  zero.  Hence 
d T 

that  at  ~ ® possible,  and  also  that  each  term  of  the  sum  (3.20) 

is  non-negative.  This  is  a very  basic  and  a non-trivial  assumption 

d T 

For  the  time  being  we  shall  try  to  avoid  it,  assuming  ^ > 0 for 

all  t e (-»,  +<»)  . 

d T 

Since  ^ > 0,  we  can  use  implicit  function  theorem  and  express 
t as  a function  of  t: 

t = 4.(t),  C-^^(t)  = C^^(t),  q“(t)  = q“(T),  S(t)  = S(x),  etc. 

d T 

It  also  follows  from  ^ > 0 that 

, a 

— — ^ — >0  (x=l,2,...,m.  (3.22) 

3q“ 

It  is  a standard  argument  (see  Velanis  that  there  is  a set 

of  constitutive  equations  of  the  form 


otherwise  the  inequality  (3.22)  can  be  violated. 

Following  this  discussion  the  following  relations  could  be 
hypothesized: 


2 


3C 


la 


= L^e^(T) 


(3.23) 


-S 


= le  = y(t)*  L30 


(3.24) 


is  one  of  the  internal  .variables  (k,A  = 1,2,3),  * is  the  con- 
volution operation: 


f * g = / f(T-T')  g(T')  dr', 

Tft 


Lj^,  L^r  L3  are  linear  differential  operators,  which  are  convolutions 
(see  Mikusinski  [33]).  These  relations  reduce  to  usual  assump- 
tions of  materials  with  memory  (in  the  material  time  parametriza- 
tion),if  the  form  of  the  operators  L^,L2»L3  are  specifically 
given  as  in  [45I  / • 

For  purposes  of  variational  formulation  we  shall  leave  the 
constitutive  relations  in  the  form  (3.23)  and  (3.24). 

It  is  apparent  that  in  the  "material  time"  parametrization 
(3.23)  is  Hooke's  law  if  is  independent  of  t,  while*<I*j^, 


r 


1 


• ( t Lj (x))  ^ • 
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is  the  identity  operator. 

Of  course  the  equations  (3.23),  (3.24)  have  to  be  consistent 
with  the  first  and  second  laws  of  thermodynamics. 


s > 0 - Lje  > 0 


(3.25) 


- ^ L,  ^ ^ 

3a^j  1 dt  dx 


(3.26) 


- ^ --i  sf  i 


(3.27) 


(i,j  not  summed  in  3.27) 

If  we  put  e = const  = 9^,  Q h const.,  = 0,  i.e.  ignore  the 

thermal  effects,  we  obtain  mechanical  laws 


K 


®o  • 4 


(C^7*  g)  . . L a) 

dOij  1 


(3.28) 


■ 

¥ 


Qq’^2  B)  < 0. 


(3.29) 


We  can  only  hypothesize  at  this  point  that  the  constitutive 
equations  must  be  of  such  form  that  the  equations  (3.23),  (3.24) 


M 


1 


i 


or  (3.26),  (3.27)  represent  a critical  point  of  some  functionals 
*1'  *2'  that 


A*q“  = f^(C^^,q®,iJ;)  = i, 
a 


A p„  - • 


A*  = 

" dT 


A*;  H2  , 


A = 


A:  Hj  » 


^0  expect  f , cp^  to  be  convolutions.  (3.15)  is  of  course 
of  the  required  form 


(-2.)  !..  = $.  , where  we  identify  =tj;, 

P 17  4^ij  * 


corresponding  to  the  variational  principle 


8(0^^)=$.  , where  B is  a linear  operator,  which  is 

4 — i -i 


mapping  space  into  space,  with  properties  of  these  spaces 


still  undetermined. 


/ I ■ • 1 


A simplified  versions  of  such  equations  of  state!  have  been 
also  suggested  in  Russian  literature 

_ fOB  where  is  the  linear  strain  component 

■^ij  ■ “^ijaB  ' 


and  where  oj. 


1 / 3 


ai,o6 


(_L_  + -1—  ) (-1^  + i(C)  + 

aB  ' 

:ir 


^ 3C“®  3C“^  3C^^  3C- 


41 


and 


= i (-2^  > 


T-  ) 4.(0)  , 


leading  to  variational  principle  of  Hu,  corresponding  to  a critical 
point  of  the  functional; 

Il(T,e,u)  . / .V.u.),dv 


- / 

2Sl 


(p^Uj  ds  - V • T^^Uj)d§, 


(see  0.  Guz  ti7l) 


For  computational  purposes  it  is  not  good  enough  to  know 
that  certain  equations  of  solid  mechanics  represent  a critical 
point  of  a functional.  Iterative  techniques  which  have  been 
eminently  successful  in  such  computations  always  had  the  additional 
information  that  this  functional  attains  a maximum  or  a minimum. 

To  m2dce  certain  that  this  occurs  we  need  the  generalization  of  the 
concept  of  second  and  higher  order  derivatives.  This  will  be 
done  in  the  next  chapter.  It  is  not  clear  at  this  point 

how  to  formulate  "a  universal  variational  principle"  which  would 
incorporate  the  equations  of  state,  (if  it  exists?),  and  we  shall 
stop  conjecturing  at  this  point  and  return  to  an  area  where  varia- 
tional principles  are  easily  established. 


Chapter  4 , 


The  second  derivative. 


4.1.  Heuristic  comments. 

Suppose  X = . . . »x^)  , ^ = (yi»y2' • * • vectors  in  r”, 

r"'  respectively,  and  that  we  have  the  functional  relation  x^^  = 

Xj^(y),  Xj  = X2(y)  ...  i.e.  a map  R*”  .*■  r”.  What 

meaning  do  we  assign  to  - — ? In  this  case  we  have  suggested 

3y 

in  Chapter  1 that  the  n x m (Jacobian)  matrix  of  partial  deriva- 
tives / 


3X. 

<5yT> 


i = 1,2,. ..,n 
j = 1,2, ... ,m 


plays  the  part  of  the  Frechet  derivative.  If  X € Hj^,  Y € ^2' 
where  » ^2  infinite  dimensional  Hilbert  (or  Banach)  spaces, 

then  the  meaning  of  the  Jacobian  matrix  has  to  be  redefined. 

However  we  have  to  recognize  that  even  in  the  finite  dimensional 
case  X € Y € R*"  the  Jacobian  matrix  is  not  an  element  of  either 

R*'  or  B™,  but  a linear  mapping  from  B*^  into  r"*.  Its  conjugate 

transpose  tr~}  is  a linear  map  from  r”^  to  r".  To  even  attempt 

3X . 

2 

^ X 

to  define  what  is  meant  by  - — ~ we  need  a few  basic  concepts  of 

3Y-^ 

functional  analysis. 


4.2.  Tensor  products. 


For  purposes  of  clarity  we  shall  carefully  distinguish  between 
elements  of  a Banach  or  Hilbert  space  H,  and  its  algebraic  dual. 
In  the  definition  of  the  algebraic  dual  we  use  the  assumption  of 
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linearity  only.  No  boundedness  tor  continuity)  properties 
are  assumed  for  the  functionals  in  H*. 

Let  E be  a Banach  space,  and  E*  its  algebraic  dual.  We  con- 
sider a Cartesian  product  of  E*  x E*,  i.e.  all  posible  sets 
of  ordered  pairs  cp  € E* , ip  € E*  (denoted  by  ) , and 

define  the  map  E -»•  R by  the  operation  4'(Xj=  <ip,X>  which  is 
defined  for  any  X C E.  We  shall  use  the  following  notation. 

Vectors  subscripted  by  indices  ...  will  denote  elements 

of  E*',  while  superscripts  will  denote  elements  of  E.  Our  algebraic 
operation  then  assumes  the  formal  representation: 


(9^  0 


(4.1) 


where  <l/'j,X^>  is  a real  number  obtained  by  applying  the  linear 
map  'I'j  € E*  to  the  vector  X^  € E.  Hence  (cp^^  0 liij)  is  a linear 
map  whose  domain  is  E and  whose  range  is  contained  in  E,  and  is 
of  dimension  one,  or  zero.  In  the  same  manner  we  define  the 
tensor  products  (cp^  © ip^);  e*  -*■  E* 

(cp  ® iji  j ) : E -►  E* 


(q)^  0 ip^):  E*  ->■  E*  by  the  algebra! 


ic  rules: 


(cp^  0 = <4-.,  X^>qj^ 


(4.2) 
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9 . (4.3) 

Again  we  see  that  the  range  of  this  map  is  of  dimension  ^ 1. 
Higher  order  tensor  products  are  defined  in  the  scune  manner.  For 
example  0 ipj  ® rij^)  € E*  x e*  ^ E*  is  defined  as  a map  from 
E X E into  E,  or  as  a map  from  E into  E x e. 


(<Pi  ® 0 nj^)  (x’'  0 y^)  = <nj^,X^>  • (4.4) 

We  assume  associative  property  of  this  algebraic  operation  in 
the  following  sense: 


: { 


((<5Pi  8 4-j  0 nj^)X^) 

= <nj^»X^>  (cp^  0 I^j)  Y^ 

= <njj»X^>  <ipj,Y^>(p^ 

= (.Pi  0 0 Hj^)  (x*^  0 Y^)  (4.5) 

We  shall  not  assume  that  such  compositions  of  mappings  are  com- 
mutative, and  indeed  such  assumption  can  not  be  consistently 
supported  in  the  general  case.  Moreover  we  try  to  pattern  our 
discussion  to  agree  with  the  usual  engineering  ideas  of  what  a 
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"tensor"  is. 

We  find  no  reason  why  the  notation  J could  not  be  used  to 
denote  a specific  tensor  product  (cp^  0 4/^).  At  this  point  we  are 
going  to  make  further  assumption  (which  generally  may  exclude  some 
important  considerations  in  mechanics,  but  seems  justified  in 
solid  mechanics).  We  are  going  to  restrict  the  spaces  E,  E*,  (E**) 
to  be  topological  duals  of  each  other.  Hence  all  linear  maps 
defined  above  are  now  assumed  to  be  continuous.  This  could  be 
labelled  "finite  energy  hypothesis".  Hence  if  our  discussion 
concerns  Hilbert  spaces,  we  can  use  Riesz  representation  theorem. 
This  will  considerably  simplify  the  definitions  of  higher  order 
derivatives  which  are  given  below. 

Let  f:  E -►  E be  a continuous  function.  We  shall  use  the 

notation  Y = f (X) , or  Y = Y (X) , X € We  define  the  Gateaux 

derivative  of  f,  (which  also  may  be  called  the  directional  deriva- 
tive of  Y with  respect  to  X in  the  direction  of  h)  as 


lim 

t->0 


Y(X+th)  - Y(X) 
t 


■MX,h)  € E. 


(4.6) 


provided  this  limit  exists.  We  claim  that  for  any  Z€  E* , <Z , 
is  a linear  functional  of  h,  i.e.  <Z,ip(X,ch)>  = c<Z,4)(X,h)>  and 
<Z  , 4»  (X,hj^+h2)  ^ = <Z  , ij' (X,hj^)  > + <Z  , t|/ (X,h2)  > . Each  follow  easily 
from  assumption  of  the  existence  of  the  limit  (4.6)  in  some  open 
region  of  E.  If  we  assume  continuous  dependence  of  4/  on  h then 
by  Milgram  - Lax  theorem  (the  bilinear  version  of  Riesz  represen- 
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tation  theorem)  c<Z,i^>  can  be  written  as  <AZ,h>,  but  in  this  case 
the  operator  A is  the  operator  of  multiplication  by  a constant 
and 

3Y 

^ '"y"  > “ c^Z^h>» 

a^h 

Again  if  we  assume  that  c depends  continuously  on  Y,  (linearity 
is  obvious)  then  using  Riesz  representation  theorem  we  conclude 
that  there  exists  a vector  cp  6 E*  such  that  c = <cp,Y>.  Hence 


aY 

, Z>  = <i|/,Z>  = <cp»Y>  <h,z> 

3Xh 


(4.7) 


Again  using  the  Lax-Milgram  theorem  we  can  rewrite  this  product 


<’^,Z>  = <Acp,Z>  <Y,Bh> 


(4.8) 


where  A,B  are  linear  operators 

B:  E ->  E* 

A:  E*  E. 

Since  this  is  true  for  abritrary  Z 6 E* , we  have  equality 


= <Y,Bh>  ACP  = <Y,^>  ACp  = <Y,c>u,  where  u = A«P. 

o X. 


Hence  the  operator  -r^  can  be  represented  as  a tensor  product 

3Xh 


= (y  ® C) / n € E,  c € E* 


Y = <Y,C>U. 


(4.9) 


Similar  argument  gives  us  a representation 


= <h,t>y 


(4.10) 


where  C is  of  the  form  c = CY,  and  C is  a linear  operator, 

* 

C;  E -»•  E*.  In  this  entire  argument  the  domain  of  the  operators 
A,B,C,  is  assumed  dense  in  a sufficiently  small  neighborhood  of 
a region  considered  for  the  respective  vectors  h,  Y (in  E,  E* 
respectively) , and  the  existence  of  adjoint  operators  allows  the 
necessary  manupulations  of  the  products  <»>.  We  can  generalize 
the  concept  of  a gradient  by  observing  that  our  assumption  of 
continuity  with  respect  to  h,  allows  to  define  grad  Y = Y^^  = 


(y  Q 5) , since 


= <5  ,h>vi. 


(Recall  the  definition  of  a gradient  of  a functional.)  In  parti- 
cular if  there  exists  a functional  V(X,x) , x € E»  X € E,  such  that 
Y = grad  V(X,x)  € E*  (or  Y = V ) , then  we  can  denote  by 


V^X  = (p  9 U = Vx- 


Si' 


g 


■ ■ 


If  we  assvune  that  V = V(X,x)^  X,x  € E,  then  a theorem  of  Valn- 
berg  states  that 


V = V . 
''Xx  ''xX 


(4.11) 


(This  does  not  imply  (w  Q C)  = ( £;  ® w)!)  (see  [49J  » chapter  Z ) 
In  fact  commutativity  of  the  "second  derivatives"  is  a necessary 
condition  for  the  existence  of  the  "potential  functional"  V. 

The  idea  of  second  derivative  is  easily  generalized  to  the 
case  when  V:  E x E*  -*•  R is  a functional  depending  continuously 


on 

..,Xj^€E  Xj^+i/Xj^+2' • * * '^n 

6 E' 

E = 

El  0 E^ 

0 . . . 9 Ej^,  E* 

9 E*. 

and 

Xi  e 

El...  X„  € 

^l'^2' * ’ 

• '^k 

^Vl---^*n 

of  E* , where 

...  V 

can 

E*. ...E*  of  E*,  where  ...  V can  be  regarded  as  genera- 

K+X  n A,  ^ 

1 n 

lications  of  the  idea  of  partial  differentials.  Consider  a 
functional  V(X^,Yj)  (not  necessarily  linear)  whose  values  in  R 
depend  on  X^  € E,  Yj  € E* . Then  provided  all  the  derivatives  shown 
exist,  we  can  construct  the  following  diagram  ; 
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Figure  4 . 1 

It  shows  the  increasingly  higher  order  tensor  products  obtained 
by  performing  repeatedly  the  gradient  operation. 

Positive  (positive  definite)  second  order  tensor  product. 

Since  (c  ® n) , c € E,  n € E,  (or  (p  0 v)  p€  E*,  v€  E*) 
can  be  regarded  as  maps  from  E*  into  E*  (or  from  E to  E)  we  can 
define  the  concept  of  positive  property  of  such  map  by  the  usual 
definition:  (^  0 n)  C/n  € E is  positive  if 

< (C  ® n)  Y,  Y> 


> 0 for  all  Y 6 E*. 


(4.12) 
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Similarly  (u  0 v)  , v (:  E* , v £ E*  is  positive  if 

< (u  0 v)  X,  X>  ^0  for  all  X e E.  0 n)  is  called  positive 

definite  if  < (^  0 n)  Y,  Y>  >0  for  all  Y£  E* , and  <(C  0 n)Y,Y> 
= 0 implies  Y = £ E* . Similar  definition  is  given  for  (p  0 v)  , 

V e E*,  V C E*. 

Vainberg's  lemma  and  some  of  its  consequences. 

Let  Ej^  be  a subspace  of  a Banach  space  E,  and  V a functional 
V:  E ->  R,  the  value  of  V depending  on  a vector  X 6 Ej^,  and  on 

Y € E2,  where  E^  0 E2  ^ E,  E^  f1  E2  = Suppose  that  V is  defined 
in  some  neighborhood  of  Xq  € Ej^.  Then  a necessary  condition  for 
Xq  to  be  an  extremal  point  of  V,  when  V is  restricted  to  E^^,  is 

= 0,  V|„l  . (I  means  restricted  to).  A sufficient  condition 
for  a min.  (max)  of  is  that  = 0,  is  positive  definite 

(negative  definite) . 

EXcimplos  of  applicatio:'! . 

Consider  the  bohavior  of  the  (non-linear)  functional 

V(X,P)  = <AX,X>  »<X,X>"^  + 4i(P)  X€  E ^ ^2  where 

are  Hilbert  spaces,  -■*  . 

Let  P be  fixed  (in  H2) • Find  sufficient  conditions  for  a 
local  minimum  of  V (in  H^^)  . We  compute  the  gradient  of  V in 

= 2<X,X>"^  • (AX  - V(X)X) 


Hence  a necessary  condition  for  an  extremum  of  V to  occur  at 


Xq  € is:  AX  - V(X)X  = 0 when  X = Xq,  that  is  AXq  = AX^, 


where  A = V(Xq),  which  means  that  Xq  is  an  eigenvector,  and  V(Xq) 


the  corresponding  eigenvalue  of  A. 


In  what  follows  let  us  assume  that  the  multiplicity  of  A, 


is  one,  and  Aq  is  the  lowest  eigenvalue  of  A.  The  second  deriva- 


tive is  computed  as  follows: 


(<n  9 tli)h  = V„^h  = lim 

t-<-0 


Vj^(X+th)  - Vj^(X) 


= 2t“^  (<X+th,  X+th>”^  (AX  + tAh  - V (X+th)  [X+th]  ) 


- <X,X>  (AX  - V(X)X)  } 


lim  2t"^  { + tAh  - V(X+th) (X+th)  _ AX-V(X)Xj 

t->0  <X,X>  + t <h,X>  + t^<h,h>  <X,X> 


Having  assumed  the  continuity  of  V,  we  obtain  as  t -*■  0 


V V,  - 2 (Ah  - V(X)h) 
^xx"  “ <X,X> 


Hence 


A - V(X)i 


V =2  

XX  <X,X> 


Hence  at  the  point  Xq  where  V^^(Xq)  =0,  we  have 


1 
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I 


A - V(X.)I 

^xx  “ 2 -T-rx-  > ■ 


where  I is  the  identity  map  ^ 
or 


''xx«=>'o>  ' - *o^>  • 


Hence  the  extremum  of  V at  Xq  occurs  on  a small  neighborhood  of 

^ c ri.Ly  if  for  every  vector  c € N„  it  is  true  that 

-^0  ^0 

f ^ ^ ^0 

<A?,C>  - An<C/C>  - <Vy(c),C>  >0  •< 

1 c /a 

Or  <Ai,i>-  yic<5,5>-<V^f,S>  *=° 

In  particular  if  A is  positive  definite,  completely  continuous, 
and  Aq  is  the  lowest  eigenvalue,  then  the  Rayleigh  quotient  V 
attains  its  minimum  value  at  Xq,  such  that 


V(Xq)  = Aq  (since  <Vj^(Xq),  Xq>  = 0). 

It  is  clear  that  in  a sufficiently  small  neighborhood  of  Xq, 
is  positive  definite  and  a local  minimum  takes  place  for  V(X)  at 
Xq.  a global  theory  is  much  harder,  and  it  is  unreasonable  to 
expect  that  the  signs  of  first  and  second  abstract  derivative  in 
some  neighborhood  should  imply  anything  about  the  global  behavior 
of  a given  functional. 
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An  Example;  Thin  plate  theory. 

The  plate  satisfies  the  linear  equation 


(D(x,y)  w(x,y))  - (1-v)  ^^(D,w)  = g(x,y)  (4.13) 


in  a region  n c R 


2 2 2 
a4,»  ov  3 a 3 B _ 3 
(Here  ^ (A,B)  = — ~ — 2 ” ^ 

3x  3y 


2 2 2 2 
3 A 3 B . 

3x3y  3x3y  3^2  ^^2 


with  boundary  conditions  w = = 0 on  3Q.  Where  the  boundary 

o ri 

3fi  of  n is  assume  smooth,  except  for  a finite  number  of  external 
corners.  (i.e.  corners  like  the  one  shown  below  are  not  permitted.) 


figure  4.2. 


This  can  be  recognized  as  the  Euler-Lagrange  equation  for  the 
functional 


I = / {D(V^w)^  - D(l-v)  (0^(w,w))}  dxdy 


(4.14) 
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It  is  not  very  convenient  (or  very  practical  in  most  engineering 
applications) 

to  assume  the  existence  of  all  derivatives  appearing  in  the  formal 

Euler-Lagrange  equations  (4.14)  for  the  functional  4>(w)  . Physically 

we  only  require  twice  differentiability  of  w,  and  L2  property  of  the 

2 

derivatives  (w  € H (n)).  We  introduce  the  following  maps 

N.;  - H^(n) 


where  N is  the  following  positive  definite  matrix 


D(x,y)  > 0 in  and  0 < v < j for  physical  reasons. 

T is  a 2 X 4 matrix  for  first  order  differential  operators^T*  is 
the  transpose  of  T 


v.j 

I 

i 


j. 

V 

it, 

r- 


It  can  be  checked  that  = N,  and  is  again  a 

1/2 

positive  definite  (and  invert: hie)  matrix,  and  fortunately  N 
is  symmetric.  The  mapping 


A;  ->■  L2(n) 

has  a formal  adjoint  A*  L2(n)  -^^*(n)»  and  Aw  =?/tt  A^^=  q 

(4.17)  where  the  derivatives  now  become  distributional  deriva- 
tives, L2(n)  being  regarded  now  as  a subset  of  5*(n),  where 

= Cq(£J)  . 

A*  = div  (T*N^/^) , (4.18) 

i . e. 


56  a 


A*itl=  div  (T*  ) = div  (T*M)  = q(x,y) 


(4.19) 


The  distributional  character  of  q(x,y)  being  supported  by  the  usual 
interpretation  of  admitting  point  loads,  or  point  moments  in  en- 
gineering practice.  Hence  we  have  the  following  mapping  diagreun 


a* 


q 


where^*^tw  = q is  the  original  differential  equation  of  static 
deflection  of  the  plate  : 


V^(DV^w) 


- (1- v)  ^ (D,w)  = q. 


(4.13) 


We  can  observe  that  following  our  theory  in  chapter  2,  the  set 


aw 

* 

= q . 


(4.20®) 

(4.20^) 


represents  two  variational  principles.  We  introduce  the  Hamil 
tonian 

W = <q,w>  + 


[A,B}  = / (A^Bj^  + A2B2  + A3B3  + A^B^)  cUdy 
0 


where 
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while 


<a,b>  = J (a  • b)dxdy. 
C2 


The  dual  variational  principles  given  for  the  static  plate 
problem  have  been  given  in  [25]  by  the  author,  and  may  be  summarized 
in  the  observation 


w„  = q 


(4.21^) 

(4.21^) 


where  denotes  the  gradient  of  W restricted  to  the  Hilbert  space 
(with  the  product  while  is  the  gradient  of  W restrict- 

ed to  H (with  the  product  <,>.  The  examination  of  second  Frechet 
w 

derivatives  i.e.  of  the  tensor  products  (W_)^  , (W  ) reveals 
that  the  corresponding  Lagrangian  functional 


L = {Orw,W}  - W (4,22) 

attains  a minimum  over  the  admissible  w € H^(n)  with  nt  regarded  as 
fixed,  and  a minimum  over  the  admissible  L2(n)  with  w regar  - 
ed  as  fixed  at  the  critical  points  corresponding  to  the  solution 
of  a system  of  equations  (4.20^),  (4.20^).  Combining  this  system 
into  the  form  AA*w  = q gives  a single  equation  (4.13)  which  is 
the  basic  static  deflection  equation  of  thin  plate  theory. 
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The  dual  variational  principles  given  here  were  originally 
discovered  by  the  author  in  [Z5]^  At  this  point  we  wish  to  make 
a more  detailed  examination  of  the  equation  (4.13)  and  find  alter- 
nate variational  formulation  for  this  equation. 

We  have  the  following  mapping  diagram: 


^1/2^  j.^1/2 


figure  4.2 


Let  us  concentrate  on  the  following  portion  of  the  diagram  (4.2) 


G is  the  Green's  function,  or  influence  function  in  engineering 
terminology;  w = G * q,  where  * is  the  convolution  integral. 

^ is  the  corresponding  map  ^ * Q = vw  = S»  which  is  to  be  determined, 
but  whose  existence  is  not  hard  to  prove,  if  we  assume  the  existence 
of  the  Green's  function  G(x,y).  (We  do  not  write  G(x,5,y,n)  since 
the  convolution  product  takes  care  of  introducing  the  "translated" 
variables  C , n . ) 
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The  dual  equations  corresponding  to  the  diagram  (4.3)  are 

N^/^ts  W (4.23®) 

Q = Wg  , (4.23^) 

except  that  W is  a functional  defined  in  different  pair  of  Hilbert 

(2) 

spaces:  which  is  the  same  as  before,  and  Hg  = ( fi)  if 

point  moments  are  excluded,  and  only  point  loads  allowed  as  loads 
which  are  not  represented  by  functions.  Note;  For  a classification 
cf  admissible  (distributional)  loads  of  thin  plate  theory  see 
author^  paper  [3/]. 

If  point  moments  are  allowed  Hg  can  not  be  embedded  in  L2 ( n) 
and  is  not  a Hilbert  space,  but  only  an  inner  product  space  which 
can  be  at  best  regarded  as  a Rigged  Hilbert  space  in  the  termino- 
logy of  Gelfand  and  Shilov  (see  [14- ] volume  4,  chapter  1,  section 


(2) 

4).  For  the  purpose  of  this  discussion  we  identify  Hg  with  L2  (fi) 
ifipeo 

with  thet  product 
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The  generalized  Hamilton's  canonical  equations  are 


T*N 


1/2 


= Q = W 


S ' 


(4.23®) 


(4.23^) 


L = I {Vl,iri}  + (S,Q)  = {N^/^TS,»l}  - W , (4.24) 


'H, 


are  exactly  the  equations  (4.23-),  (4.23-)  expressed  as  restrictions 
of  appropriate  Frechet  derivatives  of  L to  the  spaces  H^and  Hg 
respectively. 

The  tensor  products  ^slfC  ^SS  arranged  in 

a 2 X 2 matrix  form  of  operators; 


T*n^/2  ^*liT 


uiffu  ms 


w 


s??t  '^SS 


(I  is  the  identity  operator.)  V^g  = (W^^^  * is  a necessary  condi- 

Si  b 

tion  for  integrability  of  the  differential  system  (4.23  ),  (4.23  ), 
while  positive  definite  nature  of  I:  T*NT:  Hg  -►  Hg 

assure  the  existence  of  a double  minimum,  i.e.  L attains  a minimum 
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in  H^(for  a fixed  S)  if  the  equation  (4.23^)  is  satisfied, 
and  a minimum  in  Hg  ^for  a fixedT?^)  if  the  equation  (4.23^) 
is  satisfied.  These  variational  principles  suggest  numerical 
techniques  of  the  type  introduced  by  Greenspan  in  [C,  \ for 
computation  of  approximate  solutions,  and  these  principles 
may  be  simpler  than  the  known  principles  suggested  by 
equations  (4.21^),  (4.21^).  The  advantage  of  obtaining  a 
symmetric  form  of  corresponding  operators,  and  of  not  having 
the  matrix  N appearing  in  only  one  set  of  equations  are 
considerable  in  actual  computation.  The  positive  definite 
nature  of  the  matrix  N allowed  us  to  restate  our  problem  in 
terms  of  rather  than  M.  Physical  limitations  on  such  repre- 
sentation are  clear:  D(x,y)  > 0 for  all  x,y  € ?f,  (we  do 

not  have  holes  in  n.  If  we  do, let  us  relabel  what  the  region 
n really  is!),  and  ^/2  > v > 0 which  is  readily  recognized  by 
any  engineer  .,s  the  expected  behavior  of  a physical  solid, 
namely  that  the  solid  does  not  expand  volumetrically  under 
pressure. 
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chapter  5 . 

Boundary  Value  Problems 

5.1  Integration  by  parts  formulas. 

The  geometric  theory  associated  with  differential  forms  is  due  to 
Grassman.  The  expository  texts  include  H.  Flanders  1^1,  P.K. 
Rashevskii  "Geometric  theory  of  partial  differential  equations 
(in  Russian) . The  basic  axiom  for  exterior  product  of  two  differ- 
entials is 


dx  . Adx  . = -dx  . Adx  . , 
1 3 3 1 


where  the  symbol  A is  frequently  omitted.  Denoting  by  — ^ = D{ 

oY 

the  Jacobian  tensor  product  corresponding  to  a coordinate  transfor- 
mation, we  can  derive  the  general  formula 


.1 


a . . dx,  dx-  . . .dx, 

1 2 K 


I « . 


. 1.  . ^ 

3i<32--<Dk 


% -i  •••^yv 

3 1 • • • Jk  1 ^ 


: <,? 


where  a ' . 


^i  =>k 


/ X . . . . X . 

^1  ^k 


\ 


! 


.. 


\ y ■ . . . y . ! *1  1C 


The  generalized  Stokes  theorem  can  now  be  stated  for  an  arbitrary 
oricntable  manifold  u,  whose  boundary  an  is  smooth. 


X >1 
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/ d a = / a (5.1) 

(See  for  example  a monograph  by  de  Rham) . The  (exterior) 

differential  forms  provide  a natural  setting  for  certain  class  of 
boundary  value  problems.  The  formula  (5.1)  embraces  all  classical 
results  of  changing  integration  over  to  a corresponding  boundary 
integral.  We  shall  not  pursue  the  subject,  and  try  to  offer 
another  exposition  of  algebras  of  differential  forms  on  manifolds. 
We  shall  only  assume  that  appropriate  formulas  for  integration  by 
parts  exist,  and  can  be  deduced  in  the  special  cases  by  the  use 
of  formula  (5.1) . 

We  first  offer  an  analogous  definition  of  a Frechet  derivative 
of  a functional  $ (X)  such  that  the  domain  of  the  functional  $ is 
the  union  of  tv;o  Hilbert  spaces  of  functions,  i.e.  X € U ^2, 
where  is  a Hilbert  space  of  functions  v^hose  domain  is  a set 
with  a product  ( , ) ( qj  » snd  H2  is  a space  of  functions 
whose  domain  is  dQ,  (the  boundary  of  n)  the  inner  product  in  H2 
deing  designated  by  ( , ) .„•  Again  the  Frechet  derivative  in  each 
space  and  H2  is  defined  by  the  formulation  identical  to  section 
1 of  chapter  1.  The  formulation  of  a gradient  of  | in  and  H2 
presents  no  problems,  provided  the  class  of  functions  which  we 
have  in  mind  restricted  to  n and  to  3Q  respectively  do  form  ap- 
propriate Hilbert  (or  at  least  Banach)  spaces.  A very  serious 
problems  of  formulation  arises  when  boundary  values  are  not 
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functions  (and  possibly  not  even  distributions  over  some  “reasonable" 
test  space) . For  the  time  being  we  shall  assume  that  this  is  not 
the  case  and  all  possible  boundary  values  are  functions  which 
belong  to  a Hilbert  space  H2(an).  Whenever  we  use  a symbol  ft 

U H2»  we  mean  that  the  same  function  f : n U an  -►  R is  considered 

simultaneously  as  an  element  of  if  x € is  a point  of  n»  and 
of  H2  if  X € ^ ^ is  a point  on  an. 

If  f is  not  sufficiently  smojth  in  the  neighborhood  of  an» 
then  the  integration  by  parts  formulas  which  form  the  special  case 
of  (5.1)  are  not  applicable,  and  the  behavior  of  products 
and  (f,cp),^  is  completely  unrelated.  Hence  certain  continuity 

0 A4 

(and  smoothness)  properties  must  be  specified  for  functions  in  Hj^ 

near  if  the  formulation  of  boundary  value  problems  is  to  have 

a unique  solution.  To  offer  a trivial  example  we  consider 

solutions  for  the  Dirichlet's  problem  in  the  unit  disc  of  the 

complex  plane,  looking  for  y L2(3fi)  solutions  with  no  continuity 

2 

requirements.  Suppose  we  wish  to  solve  the  problem  v u e 0 in 
(the  open  unit  disc)  u = 0 on  an  (the  unit  circle) . In  the  class 
of  differentiable  solutions  continuously  approaching  the  boundary 
value  there  is  only  one  solution  u = 0 in  n.  If  no  continuity  is 
required  on  approaching  3n,  then  u e 5j  in  n , u = 0 on  an  is  another 
possible  solution  among  infinitely  many  candidates  for  the  solution 
of  this  problem.  In  fact  it  is  analytic  in  n.  However  the  possi- 
bility of  discontinuity  on  approaching  the  boundary  made  a speci- 
fication of  the  boundary  condition  completely  worthless.  This 


fairly  trivial  example  interprets  the  meaning  of  our  previous 
remark  that  in  general  the  behavior  on  SI,  and  the  behavior  on 
of  functions  in  the  class  L2 ( fi)  U L2(3fi)  is  completely  unrelated. 

We  have  to  digress  into  discussion  of  the  following  basic  question. 
Let  B2  be  some  Banach  spaces  of  functions  whose  domain  is 
n ^ r”.  Let  A be  a given  operator  A:  Bj^  B2.  Let  f € B2  be 

given.  Let  B^  be  some  other  Banach  space  of  functions  whose  domain 
is  an.  Let  <P  € B^  be  given.  Does  there  exist  a vector  x € Bj^, 
such  that  x satisfies  the  equation  Ax  = f € B2  and  that  the  values 
of  function  x:  n->-R  converge  (in  some  sense)  to  <p : 3n-*-R,  in  some 
neighborhood  of  an  in  n.  That  is  can  f be  continuously  extended 
to  n U an  so  that  (in  some  sense)  it  coincides  with  cp  on  an?  In 
the  problems  of  mathematical  physics  "in  some  sense"  usually  is 
understood  to  be  pointwise  convergence  almost  everywhere.  Usually 
the  physical  interpretation  requires  that  a solution  of  such  problem 
should  exist  (since  some  physical  process  is  going  on,  representing 
a solution) , and  frequently  we  would  like  to  assert  that  such  a 
solution  should  kc  unique.  If  a unique  solution  does  exist  to  this 
problem  it  is  called  well  posed,  or  properly  posed  boundary  value 
problem.  An  entirely  different  question  is  the  continuous  depen- 
dence (stability)  of  such  solution  on  the  boundary  data  (function) 

9P  : an  -♦  R.  The  common  question  is  will "smal  1* variations  in  cp 

result  in  "small"  variations  in  the  solution?  Of  course  mathe- 
matically the  above  sentence  makes  little  sense,  but  it  can  be 
easily  translated  into  a rigorous  statement  (given  e >0,  there 
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exists  6 > 0 such  that  5 - neighborhood  of  Cp  in  is  mapped 
into  e ~ neighborhood  of  x in  under  the  inverse  map  (f,(p)  -+•  x, 
with  the  existence  of  such  inverse  map  assured  by  well  posedeness 
of  the  problem.)  The  well  posedeness  of  the  boundary  value  problem 
implied  the  existence  of  an  inverse  operator of  the  operator 
15  ! X -*■  (f»V)  € B2  U i.e.  ofj^  (f,cp).>  x.  The  stability 

in  the  sense  given  above  implies  that  fb~^  is  a bounded  operator. 
Boundary  value  problem  which  is  well  posed  and  stable  is  called 
in  the  literature  properly  posed  in  the  sense  of  Hadamard.  For 
a classical  example  of  a problem  which  is  unstable  and  therefore 
improperly  posed  in  the  sense  of  Hadamard  consider  the  Laplace 
equation 

U . 0 

3x“^  ay 

2 

in  half  plane  y ^ 0 R )»  with  boundary  conditions  ^jy_Q  = cp  (x) 

(x)  , -00  < X < +00 , where  cp  c C(-:o, +<*>).  (See 

y y=0 

Mikhlin  [32],  chapter  9,  § 5 for  a discussion  of  this  problem, 
and  chapter  25  for  a general  discussion  of  the  well  posedness 
of  boundary  value  problems  for  partial  differential  operators 
of  mathematical  physics.)  The  main  thrust  of  this  discussion 
is  the  following  conclusion:  the  boundary  value  problem  is 

not  well  posed  or  ill  posed  "by  itself".  The  well  posedeness 
depends  on  the  topologies  of  the  spaces  B^^ , B2 , B^,  and  the 
problem  may  be  well  posed  for  some  choices  of  Bj^ , B2,  B^,  and 
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^ ill  posed  for  other  choices.  In  particular  we  must  answer  the 

I problem  of  existence  and  boundedness  of  the  operator^ 

I B2  U From  this  point  of  view  it  is  easy  to  interpret 

; the  following  statement  which  physically  sounds  ridiculous,  but 

is  mathematically  correct.  For  certain  class  of  boundary  func- 
tions (analytic  of  certain  slow  growth  rate)  the  backward  heat 
conduction  problem  is  well  posed.  Here  the  boundary  is  the  phy- 
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sical  space  boundary  of  the  solid  body  and  the  time  hypersurface 
t = tj^  in  Hence  for  certain  choices  of  spaces  Bj^,  B2/ 

and  for  certain  classes  of  operators  A:  Bj^-»  B2  • B^  the  for- 
mulas of  the  type  (5.1)  make  sense.  Assuming  that  B2,  B^ 

are  Hilbert  spaces  and  the  domain  of  A is  danse  i”  and  there- 
fore A*:  ®2  “*■  ®1  uniquely  defined,  we  may  have  a unique 

representation  of  the  form  (u,A*v)  ^ - <Au, v > ^ = (Bu,v)  , for 

u € s ^ Bj^,  where  S is  some  class  of  functions.  We  have  defined 

the  inner  products  ( , ) , ^ in  B,  < , % , . in  B. , and  ( , ) 

tn)  1 V n)  2 3JJ 

in  Bj.  In  such  cases  A restricted  to  $ is  called  the  formal 
adjoint  of  A*.  In  the  context  of  this  definition  for  example 
the  operators  and  - are  formal  ad joints  of  each  other 


dt 


on  C^[0,11^  B,  where  C^[0,1]  is  imbedded  in  the  Sobolev  space 
H^[0,1].  We  observe  that  in  that  case  (Bu,v)  = '^'^lx=l  ~^'^lx=0* 
There  is  no  way  of  defining  uniquely  (Bu,v) if  we  attempt  to 
embed  the  problem  in  L2[0,l].  This  remark  illustrates  the  fact 
that  a setting  for  many  problems  of  mathematical  physics  is  pro- 
vided by  Sobolev  spaces.  See[|B]  for  a comprehensive  study  of  their 
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properties,  or  see  Mikhlin  [3<^]  for  an  expository  account. 

5.2  Critical  points  of  functionals  associated  with  boundary  value 
problems. 

Suppose  A is  a linear  map  A:  H(a)  -*■  where  H(fi)  is  4 Hilbert 

space  of  functions  whose  domain  is  n,  a Hilbert  space  of 

functions  whose  domain  is  the  domain  of  A is  H(f2)), 

which  is  dense  in  H(n),  and  f is  some  function  f:  H(n)  H(fi)  . 

We  formulate  the  following  class  of  functionals  on  H(n)  9 H2(3n) 

♦ (u,q>);  H(jj)  9 H2(3n)  ->  R, 

<t(u,v)  = <Au,v>^^j  - <f(u),v>j^j  + r(u,v).  1 

fi 

i 

j 

(We  have  assumed  for  convenience  that  H(r.)  H (3^2)  are  Hilbert 

2 

spaces  but  this  can  be  generalized.)  r(u,v)  is  a continuous 
functional  H_  x H_  F,  i.e.  (u,v)  : r(u)  , v|  ) -*•  R,  is  a 
continuous  functional  whose  domain  are  ordered  pairs  of  functions 
u,v  restricted  tc  3n,  and  regarded  as  elements  of  H2(3j2).  As 
before  the  inner  products  <Au,v>^and  <f,v>^^j  are  inner  products 
in  H(ji).  f is  regarded  as  a fixed  element  of  H(n),  u,  v are 
regarded  as  independent  of  each  other.  In  a more  general  case 
f is  a map  H(n)  -►  H(jj)  dependent  of  u.  Computing  the  gradient 
of  $ in  H(fi)  9 H2(3n)  we  have 

V ' » fv Hjl  3!!)  « <«2  » «>*• 


Here  Au  - f(u)  is  a vector  in  H(Q)  and  r^(u,v)  is  a vector  in 
H2(9n)«  (Since  we  have  assumed  a Hilbert  space  structure  for 
H(q)  and  H2(3fi)  there  is  no  mathematical  necessity  for  keeping 
track  of  duality.  Physically  it  is  of  course  important  to  keep 
the  gradient  in  the  dual  space  to  H(J^)  © H2(3J5),  distinguishing 
which  are  the  generalized  displacements,  and  which  are  the  gen- 
eralized forces,  and  keeping  the  physical  dimensions  correct. 

But  vanishing  of  the  gradient  implies  that  each  vector  in  the 
direct  sum  of  spaces  H(Q)  © H2(3fi)  must  vanish  separately.  Hence 
is  a simplified  notation  for  saying  that 


Au  - f(u)  = 0 in  H(n). 


(5.2) 
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with  a critical  point  of  the  functional 

<Au,v>^  - <f(u),v>^  + = 4,(u,v). 

Until  now  v is  independent  of  u,  but  quite  arbitrary.  We  can 
now  postulate  that  v is  a solution  of  the  adjoint  system 


V - 0 

3U 

in 

(5.6) 

iaiHl  . V = 0 

3U 

on  an. 

(5.7) 

i 


That  is  V is  chosen  such  that  ij>^  = 0,and  that  the  given  boundary 
value  problem  corresponds  to  a dual  critical  points  of  $ , provided 
A*  is  a (true)  adjoint  of  A.  This  means  that  g(u)  is  so  chosen 


that  for  all  u,  v € S,  <g(u),v>^jj  = 0,  and  <Au,v>  ^ - 

2 

<u,A*v>  (H(jj))  = 0. 

The  problem  which  frequently  arises  at  this  point  is  that 


A*  (the  adjoint  of  A)  may  not  be  uniquely  defined  in  H(n)  ® H2(3fi)f 
and  additional  information  may  be  needed  concerning  properties 


of  the  (physically)  admissible  elements  of  H(n)  in  the  neighbor- 
hood of  Our  previous  remarks  concerning  formally  adjoint 

operators  A, A*.  (H(fi)  ->■  H(n))  are  applicable  heie.  Presume  that 
h,  and  A*  are  formal  ad joints  mapping  a subset  S of  H(n)  into 
H(n),  such  that  for  any  u ^ S c H(£J).we  have  an  appropriate 


formula  (in  H ® H2) 
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<Au,v>^jjj  - <u,A*v>^jjj  = B(u,v)gjj.  (5.8) 

r 

We  see  that  the  existence  of  the  appropaate  functional  , 

().(u,v);  (H(n)  ® H2(3n))  ® H2(3fi))  -►  R 

whose  critical  points  are  represented  by  (5.4),  (5.5)  and  by 
some  equivalent  conditions  of  the  form  (5.6),  (5.7)  is  now  harder 
to  determine.  (When  A*  was  a true  adjoint  of  A we  wrote  down  the 
appropriate  functional  4>(u,v)  without  any  trouble?) 

We  have  the  relationships: 

Au  = f(u)  in  £2, 

g(u)  =0  on  3£2, 

or  more  generally  g(u,v)  = 0 on  3£2.  And  for  all  v € H 

<Au,v>(^)  - <u,A*v>j^)  + B(u,v) (3^)  = 0 (5.9) 

where  subscripts  (n),  (3 £2)  denote  the  inner  products  in  H(£2)  and 
H2(3£2)  respectively.  Denoting  as  the  Hamiltonian  the  functional 
W(u,v)  : 


<f (u) ,v>^  = W(u,v) 
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we  have  the  relation 

= f(u) 


while 


in  R , 


W = IliHiv  , where  the  meaning  of  has  to  be  defined 
u 3u  ' ^ 3u 

analogously  to  section  4. 

P1  - 

The  corresponding  Lagrangian  functional  is  denoted  L and  is 
defined  by  the  relation 

L(u,v)  * <Au,v>  - W 

* <u,A*v>^jjj  - W + B(u,v)^gjjj  . 

Here  the  same  symbol  v is  used  to  denote  v as  an  element  of 
^(S2)  ^2(3n)  • 

We  seek,  however,  a "true"  Lagrangian  L(u,v)  : H ® H2-*  R» 

given  by  a relation: 

L = L(u,v)  + C(u,v) 

C(u,v),„=  X Hj  B, 
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i.e.  C(u,v)gjj  is  a continuous  map;H2(3fi)  X H2<  3fi)  -*•  R. 

The  "true"  Lagragian  L(u,v)  should  correspond  to  a multiple 
variational  principle.  In  this  case  since  only  two  independent 
(vector)  variables  are  present  we  should  have  the  dual  critical 
points: 


Au  = f (u) , g(u,v)  = 0 can  be  replaced  by  the  following  equa- 
tion in  H and  H2  respectively. 

Au  -W^  = Ly  €H  (u,v|j^)  (5.10a) 


g(u,v)  = €^2  (5.10b) 

This  suggests  that  L(u,v)  should  have  the  following  form. 

L=  ^u,v>j,  - W(U,V),J„ 

=v|3Sl“  (5.11) 

'■u  = **''(0)  »<“'''>u(3!))  - "u(a)  * '=uO0) 

Hence  ^ (in  H « H2)  expresses  the  initially  given  equa- 

Au  * f (u)  and  the  boundary  condition 


tion 
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The  "symmetry"  formula  (5.14)  is  written  below  in  full: 


<u,A*v>  -<«rW^^v>^  + <u,C^^v>jj^ 


= <Au,v>jj  - <W^^u,v>jj  + <B^^u,v>3„  + . 

(5.14a) 

The  formula  (5.14)  was  suggested  by  A.  M.  Arthurs  in  [2  ], 
who  also  has  given  the  following  set  of  sufficiency  conditions 
for  the  existence  of  the  Lagrangian  functional: 


(5.15a) 


‘“'V'OS!)  = "=va"''"(5!l) 


(1.15c) 


The  condition  (5.15a)  is  a necessary  condition  for  the  exist- 
ence of  a Hamiltonian  function  W(u,v)  whose  domain  is  H(5i). 
(Not  H 9 H2 : ) . 

The  condition  (1.15b)  implies  the  following  equality: 
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= g(u,v)  = d on  an . 

On  the  other  hand  the  requirement 

(in  H e H2) 

requires  the  vanishing  of  vectors 

A*v  - in  H(n)  , 


and  of 


(5.10c) 


(5.13a) 


B(u,v)^ 


(5.13b) 


The  "independence  of  path"  condition  of  Vainberg  for  the  exist- 
ence of  L requires  that  the  tensor  products  L and  L are 

uv  vu 

adjoints  of  each  other,  that  is 


(5.14) 


for  any  admissible  pair  (u,v)  u,v  ^ S - H ® Hj,  implying  that 
u,v  satisfy  the  required  continuity  conditions  on  30,  such 
that  the  formula  (5.9)  is  valid. 


However  the  "symmetry"  condition  is  satisfied 

if  admissible  u,v  € S c h (with  suitable  continuity 
conditions  in  some  neighborhood  of  3q) 


<u,A*v> 


This  implies  that  C(u,v) 


must  satisfy  the  relations 


The  existence  of  the  "boundary  potential"  C(u,v)  such  that 
^vu  “ C^v  in^plies  (g(u,v))„  = 
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Since  H2(an)  is  a Hilbert  space  C^,  g(u,v)^, 

are  all  linear  maps  H2  H2,  and  the  above  equalities  make 

sense. 

Hence  under  these  assumptions  the  secessary  conditions  for 
the  existence  of  the  Lagrangian  L:  H ® H2  -»•  F are: 


(1.16) 

(1.17) 

(1.18) 
(1.19) 


w 


uv 


(fi)  = W. 


vu 


(£J) 


B(u,v)  I an  =<Bvu'^'^^an 


g(u,v)  = (an) 


c 


u 


<u,Bu^v> 


an' 


If  C(u,v)  , (?::)  can  be  chosen  to  satisfy  all  conditions  (1.16) 
-(1.19)  then  of  course  we  have  nothing  to  worry  about. 

In  general  however  this  is  not  a convenient  set  of  conditions, 
and  it  can  be  readily  checked  that  choosing  for  example 
C = -3,  we  end  with  the  original  condition  suggested  by  the 
author  in  [74]:  g i '•  P'^f’PibJp  f-o 

satisfy  only  in  special  cases.  To  restate  x eal isLloa 1 ly , tin 
problem  of  finding  a functional  L;  (H  ® K2)  ->■  R for  quite 
arbitrary  boundary  condition  g(u,v)  ==  0,  we  have  to  sacrifice 
some  conditions  on  the  list  (1.16)  - (1.19). 
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The  problem  can  be  restated  as  follows; 


Does  there  exist  v € H © H2,  and  C:  H2  -»■  R such  that 


A*v  ” in  H ( n) 


Cu  = 


Cy  = g(u,v) 


in  H2  ( 3n) 


where  W;  H ->■  R is  given  by  W = <f(u),v>^^j  ? 


The  answer  to  this  question  depends  on  the  nature  of  the 
boundary  conditions  and  on  the  operator  A. 


If  such  V,  C(u,v)  can  be  found, then  the  problem  of  finding 
a solution  to  the  given  boundary  value  problem  is  replaced  by 
a pair  of  variational  principles  for  the  corresponding 
Lagrangian  functional  L(u,v);  H © H2  ->  R. 

We  observe  at  this  point  that  if  L(u,v)  can  be  found,  then 
the  tensor  products  1*^^,  satisfy  the  equalities 


L = (-W  + A*)  + (B  + C ) 

uv  ' uv  'n  uv  uv' an 


L = (-W  + A) ^ + (C  ) , 

vu  ' vu  ' n vu  an 
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L = (-W  )„  + {(C  + B)  ), 

uu  uu  uu  an 


^vv  ” ^^vv^an 


The  remarks  in  the  previous  chapter  explain  the  following 

necessary  conditions  for  an  extremal  behavior  of  L(u,v)  at 

the  point  u = u , v = v . 

0 0 


a)  is  the  adjoint  of 

b)  The  signs  of  the  operators  are  defined 

and  remain  constant  in  some  neighborhood  of  the  point 
Vq,  v^  € (K  e H2)  X (H  ® H2). 

Examples  of  applications  of  multiple  variational  principles, 
and  of  boundary  value  problems  posed  in  this  manner  will  be 
given  in  the  next  chapter.  Many  examples  of  dual  variational 
principles  can  be  found  in  the  monograph  of  Arthurs  [i  ] and  in 
the  papers  of  the  British  applied  mathematicians  who  have  pur- 
sued this  topic  using  a largely  heuristic  approach.  See  for 
example  Arthurs  [2  ],  or  80 j. it  I ] ^j-or  L 

1 ^ 


V. 


nf  • hf’  rticiril  ,lifrifii1l  |>rri|,  | i ,n;;  f'lii' •»>  r 11  i ii' I '•  1 n r; « i f i i»in 

of  the  operators  or  of  the  boundary  conditions  are  still  largely 


unsolved. 


-1- 


Chapter  6 . 


U-i  ■ ' 

I iu 


Some  applications  of  necessary  conditions  for  the  existence  of 
a critical  point  of  a functional  to  the  theory  of  ordinary  differ- 
ential equations . Many  problems  of  physics  are  mathematically 
modelled  by  differential  equations,  which  under  soane  assumptions 
may  reduce  to  a system  of  ordinary  differential  equations  of  the 
form  jL  (W  (x)  ) + f (x  t)  - «f(t)  (6.1) 

where  L is  an  ordinary  differential  operator  which  need  not  be 
linear,  w and  f are  vectors  in  a Hilbert  space  (of  functions) H. 

The  domain  of  functions  in  H is  n-dimensional  Euclidean  space  R*'. 

As  an  example  we  shall  consider  a one-dimensional  case  studied 
by  Poincari. 


(a(t)*'^(t))'  = c(t)*f(x)-  Jp(t) 
'p(x)  -x'  = - '^(t) 


(6.2) 


Physically  the  system  of  two  equations  (6.2)  interprets  the 
motion  of  a mass,  which  may  decay  with  time  in  a force  field  which 
is  also  time  dependent. 

2 

Interpreting  p as  the  angular  momentum  p = p^=  a(t)"r  *0  , and 
introducing  a constraint  r = f ( e ) , we  obtain 


p = a (t)  [ 4'  (0  ) • e ' 

9 

(Pq  )'  = -c(t)  ?(r,e  ) = - c(t)  f(  e) 


K (6.3) 


which  is  identical  with  the  system  (6.2). 

Combining  equations  (6.3)  into  a single 


(t)  t 'f'  ( 0)1 


lations 
'}  '* 


c (t)  f ( 0 ) = f (t)  , 


an  equation  of  the  form  (6.1). 


equation,  we  obtain 
t^to, (6.4)  which  is 


For  physical  reason  it  is  usually  assumed  that  a(t)  >0 


t F.  (t  . c(t)  may  assume  negative  values.  We  also  assume  that 
0 

(0)  ^ 0. 

We  recall  that  in  the  nonlinear  motion  studied  by  Poincare 
the  physical  meaning  of  ij)(0)  was  the  distance  from  the  origin 
Moreover  in  this  case,  we  also  observe  that  p(0)  is  undefined  at 
ij<(0)  = 0,  and  the  motion  must  not  pass  through  the  origin.  Hence 
we  assume  that  ii»(0)  >0  and  along  the  entire  trajectory  ['*'(6)]  ^ 
is  a continuous  function  of  0.  We  can  now  redefine  p=  /a  P,  i.e. 
p(t)  = - /a  (t)  ii(x(t))  x'(t),  and  rewrite  the  system  (6.3)  as 

x‘ (t)  = [a(t)]’^2[«Mx(t))]"^  p(t) 

^/a  (t)  p(t)^  = c(t)  f(x(t))  - vf(t) 


\ 

> (6.4) 


Here  denotes  the  positive  square  root! 
We  introduce  a potential  function 


K(x) 


(6.5) 


X 

where  X„  can  be  taken  as  the  origin  if  / <()(?)  d^  exists,  or  we 

X 0 

can  put  Xq  = + ® if  / ^ 'p(C)  d^  exists.  Clearly  dK/dx  0 and 
K ^ is  defined. 

In  most  cases  the  choice  of  "the  zero  potential  level  Xq  is 
immaterial,  as  long  as  the  equations  of  motion  (6.2)  can  be 

rewritten  as  a generalized  Hamiltonian  system 
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^K(x)  = p,  (6.6-) 

a*  P = c(t)  f(x)  - ^(t),  (6.6-) 

where  • = (/a  ^ (6.7-) 

= -It-  (/i  .).  (6.7^) 


The  operators  <7  and  <7  * are  formal  ad  joints  of  each  other. 

Hence, we  have  a system  of  generalized  Hamilton's  canonical  equations 
already  discussed  in  chapter  2. 

Suppose  that  we  seeJc  wealc  solutions  of  equations  (6.o— ), 

(6.6—)  in  the  Sobolfiv  space  lot,&)  where  [o(,e)  is  some  interval 
of  time  (See  f44l , Chapter  2).  We  allow  6 = +~  in  some  considerations. 
For  obvious  physical  reasons  we  cannot  stipulate  twice  continuous 
differentiability  of  solutions,  restricting  ourselves  to^contin- 
uous  applied  force  ^(t),  since  in  many  instances  we  may  consider 
applied  (outside)  forces  which  are  only  piecewise  continuous. 

In  this  physical  setting,  we  stipulate 

p(t)  e L2  [<^/6),  9(t)  and  c(t)  f(x(t)) 

piecewise  continuous  in  [o(,B),  and  we  seek  solutions 

X(t)  e W2  [«<,3)  . 

The  finite  work  condition  is 

|c(t)  f(x(t))  'x(t)  I dt  ' - 

ot  ' ' 


Because  of  piecewise  continuity  of  f (x)  (and  of  x(t))  it  suffices 
if  c(t)  € L2  - for  a finite  interval  [(XjBl^for  the  finite 

work  condition  to  be  fulfilled. 


We  also  assume  for  a finite  interval  [-<,3] 


\ 


a (t)  > 0 for  all  t g [<^,3], 


(6.9) 


Since  It  IS  easy  to  see  how  to  generalize  such  results  to 
the  more  general  case  of  (6.11),  Here  w,,  w. , p,,  p.  are  real 

• P ^ D 

numbers.  We  remark  that  the  results  obtained  below  may  be  also 
derived  under  somewhat  weaker  hypothesis. 

Since  i^(x)  = — 0,  we  can  introduce  an  invertible 


transformation  K: 


x(t)  K (x(t))  = R (t) 


x(t)  = K (R(t)) 


(6.13) 


The  problem  can  now  be  posed  in  the  L2[''^/6]  setting  with  the  usual 
product  <u,v  > - /^  u(t>  v(t)dt.  We  also  introduce  the  inner 

product  called  energy  product: 


[u,v]  = a(t)  u(t)  v(t)  dt 


(6.14-) 


The  Hilbert  space  obtained  by  closing  this  inner  product 
space  with  respect  to  the  energy  norm  will  be  called  the  energy 
space.  II  > II  will  denote  the  L2[«<fB]  norm:  ||'^lP  ^ <u,u>,  and 

|l(  l(lthe  energy  norm:  |llu|lP  = [u,u]  . 

Ignoring  the  boundary  conditions , the  Hamiltonian  for  our 
problem  can  be  given  by 


W(x,p)  = i + c(t)  4>(K)  - <>(t)K 


-6- 


where  4>(K)  = /J  f k"^  (C)  dS,  (6.16) 

while  K(x)  = i(.(0  dK , fK"^  (K)  ^ f (x)  . 

oc 

(Note  that  defining  K(x)  = i|)(C)  dC,  or  iHUdC  would  introduce 

no  real  difficulties  into  all  subsequent  argrments) . 

The  corresponding  Lagraogian  functional  or  the  action  functional 
is  given  by 

L = W - <(JK,p>.  (6.17) 

With  the  introduction  of  this  functional  we  have  the  canonical 
system: 

= Wp  (6.18-) 

a*P  = . (6.18^) 

(It  and  d*  are  only  formal  ad  joints  of  each  other.  To  malce  certain 
that  variational  principles  of  chapter  4 and  5 can  be  directly 
applied  to  the  system  (6.18—),  (6.18—),  we  need  to  either  include 
the  appropriate  boundary  terms  (see  conditions  (6.12—)  - (6.12—)), 
or  to  checlc  that  the  boundary  conditions  are  natural.  In  this  case 
the  check  is  performed  by  integration  by  parts.  By  definition; 

<(3LK,p>  = /^I^TTET  • (g^  K(x))  • /a  ^(x)  • x'(t)Jdt 


/®ta  (t)  ,i/(x)  X*  ^ (K(x))]dt 
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[• 


If 

1% 


= a 


(t)  <Mx)  x'  K(x)|f  - ll  {K(x)  (a(t)  iMx)  x')}  dt 


= a(t)  (x)  K(x)  X I + <K,^*p> 


Imposing  the  boundary  conditions  (6.12—)  we  see  that  <2  and  a* 
are  true  adjoints  if 


a(t)  i|;(x)  x*(t)  K(x(t))  ^ = 0 

oc 


Since  t(»(x(t)),  and  a(t)  > 0 and  t ^ t^,  vanishing  on  the 

boundary  of  x'(t),  therefore  of  p(t),  or  of  K(x(t))  implies  that 
^ * is  a true  adjoint  of  2 and  the  corresponding  boundary  conditions 
at  t = <^  and  t = 6 will  be  referred  to  as  the  natural  boundary 
conditions. 

We  check  the  sign  of  the  appropriate  tensor  products; 


^KK  " ^KK  ^ ^ 


(6.19^) 


where 


A *1  ^ ^ 

f ' K)  = ^ f (x(K))  . 


3K 


^PP  ^ '^PP  ^ 


(6.19-) 


(I  is  the  identity  operator.) 


Kp 


(6.19-) 


(6.19-) 


Hence,  if  c(t)  does  not  change  its  sign  on  [<<,3]  and  provided 


that  f (K)  0 in  some  neighborhood,  a critical  point  Kq,  we  are 


assured  of  the  existence  of  a variational  principle.  The  extremal 
behavior  of  i-  corresponds  to  a solution  of  the  differential  equation 


1 

J 
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r 


(6.2).  Depending  on  the  sign  of  f'(K)  an  algorithm  can 
be  constructed  in  the  spirit  of  Greenspan's  techniques  in  [i6]  . 
Morever,  the  existence  of  such  a critical  point  of  L implies  the 
correctness  of  equations  (6.18—),  (6.18—)  which  are  equivalent 
to  the  original  system  of  equations  (6.2).  However,  the 
sufficiency  condition  for  existence  of  a corresponding  variational 
solution  of  the  problem  posed  by  (5.2)  is:  c(t),  a(t)  are  of 
constant  sign  on  (»,B^,  while  f (K)  is  Fr^chet  differentiable  in 
some  region  containing  the  possible  critical  point  K^. 

K(x(t))  6 V?2  [««»  g]  » hence  p t [»,g].  A corresponding 

function  x(t)  which  results  in  a stationary  behavior  of  L is 
t defined. the  weak  solution  of  the  differential  equation  (6.2) . 

COMMENT . 

We  note  that  the  L2  setting  of  the  variational  problem  was  not 
only  mathematically  convenient , but  also  reflected  some  very  strong 
physical  assumptions.  (It  was  a sufficiency  condition  for  finite 
work ! ) 

t V 

f.  ?■ 

[ 
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6.2  The  essential  boundary  conditions 


We  shall  consider  the  system  (6.2)  with  boundary  conditions 


W(*i)  = W., 


(6.20-) 


W (B)  = Wg, 


(6.20-) 


P(<^)  = 


(6.21^) 


P(3)  = Pg 


(6.21-) 


where  we  assume  that  the  boundary  conditions  are  not  natural,  and 
that  Cl  and  are  only  formal  ad  joints  of  each  other. 

As  before  we  formulate  the  problems  in  (real)  spaces  [•^»3] 
and  L2  [<*»B].  We  modify  the  definition  of  the  inner  product  by 
adjoining  a discrete  product:  [u,v]  = /aT^  u(^)  v(**). 

We  extend  the  domain  of  the  operator  (7  in  the  following  manner. 
For  every  f ^ K/B],  we  define 


61  : f (t) f (t) 


V t € (<t,  B) 


X f(a)  + (x-1)  f 


(X-1)  f(B)  + X f, 


(6.22) 


at  {t  = U {t  = 6) 


4r  • g(t))  t e (-,^) 
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g (■* ) 
g(3) 


{ A - l)  g(d)  + X g 


(6.23) 


g(e)  + ( A " 1)  g 


at  (t  = 0^)  U (t  = 6 }• 


f g_,-  g are  some  a priori  given  real  numbers. 

““63 

We  check  the  properties  of  the  product 

^3'f,  g^  using  integration,  by  parts; 


g}  = <af,  g>  + a [£(«<),  g(ot)]  + 


(A  - 1)  [f  g(<*)]  + (A  - 1)  [f(3),  g(3)]  + A[f^,  g(6)]  (6.24) 

P 


(f/  :7*g}  = <f,  (7*g>  + (a  - l)  [f(c),  g(a)  ] + 


Mf(-'),  gg]  + z^[f(3),  g(3)]  + (A  - 1)  [f(s),  g^] 


This  is  manipulated  in  to  the  form 


(6.25) 


{'Tf,  g}  = {f,i?*g}  + A^f(ot),  g^]  - [f^,  g(3)]j  + 

(A  - 1)  ^[f(6),  g^]  - [fg,  g(<^)]j 


(6.26) 


We  can  eliminate  the  term  with  (A  - 1)  coefficient  by  putting 
X = 1,  or  the  term  with  A by  putting  X = 0. 

Suppose  the  boundary  conditions  are  of  the  form; 


f(<)  = f 


g (e)  = g 


(*) 


3 


r 


-11- 


Then  substitution  of  ^ = 1 gives  us 


0t,  g}  = {f,(2*g}  + ^[f (•<•),  gg]  - [f^,  g(B)]j 


and  we  have  the  equality 

Cyf,  g}  = {f/i7*g} 

if  the  boundary  conditions  (*)  are  fulfilled.  Similar 
arg  ment  works  for  boundary  conditions 

g(‘)  = g^ 
f(g)  = f, 


'6 


I (**) 


where  we  need  to  put  A = 0 to  arrive  at  the  same  conclusion. 
In  the  case  of  boundary  conditions 


j.  , 

V 
c . 

T'- 


f(^)  = f. 

f(B)  = f 


B 


> (***) 


we  need  to  use  some  value  of  A,  such  that  A / 0,  A / 1. 

Since  g(<<)  and  g(B)  are  not  given, the  constants  g^,  g^  are 
also  not  given  a priori  and  we  have  to  regard  them  as  unknown 
functions.  Hence,  identifying  g with  the  variable  g(*^)  and 
q„  with  g(B)  , we  have 

P 

{frj*g}  = g}  + {A[(fK)  - f ),  g(B)]  + 


(1  - A)  l(f(B)  - fg),  gK)]},  A 0,  A 1, 
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Hence  {<7f,  g}  = {f»£7*g}  if  the  values  assumed  by  f(t)  at  t = «t. 
and  t = B are  exactly  the  values  specified  f and  f respectively. 

a,  g 

An  identical  argument  works  for  the  case  of  boundary  conditions 

g(<^)  = g^ 

g(e)  = 

We  are  ready  to  introduce  the  Hcimiltonian  and  Lagraiagian  functionals 
for  the  equations  (6.2)  with  boundary  conditions  of  either  one  of 
the  types  (*) , (**),(***),  or  (****). 

The  Hamiltonian  is  given  by 

W(x,  p)  = ^ <P»P>  + c (t)  $ (K)  - ^ (t)  K 

which  is  identical  with  (6.15).  while  the  Lagrattgian  is: 

W(x,  p)  - {K,<3*p}  = L(x,  p)  (6.27) 

If  boundary  conditions  are  of  the  form  x(o>r)  = x(q)  = 

can  replace  them  by  K(x(rf))  = K^,  K(x(b))  = since  K is  a known 

function  of  x. 

The  vanishing  of  Frechet  derivatives  of  the  Lagrangian 
functional  corresponds  to  the  following  sets  of  equations. 


W - 
P 


= 0 in  («f,6) 


(6.28) 


on  { » U { e ) 


(K(^)  - K J p(g)  = 0 


(K(e)  - K ) p(<)  = 0 

p 


Wr  - ^*P  =0  in  (*,g) 

As  before  we  check  the  signs  of  the  tensor  products 


(6.28-) 


(6.29) 


^KK  ~ f (K)  in  (^/g) 


(6.30) 


Sp  ' ”pp  " ^ 


and  we  check  the  relation: 


- 

l"  “1 

K(B)  - , 0 

p 

' A 

on 

0 K(o^)  - K 

•L 

e ! 

. 

_ J 

(6.31) 


)* 


We  can  now  state  the  following  variational  principles: 

Choosing  (arbitrarily)  x(t)  such  that  K(x(t))  satisfies  the  given 
boundary  conditions  K(x(«;)  ^ K^,  K(x  ( p) ) - v , ♦ hnn  rnn,p"» 
p(t)  to  satisfy  the  relation  =:7K,  the  actual  solution  of  (6.2) 
and  consequently  of  (6.29)  will  minimize  the  Lagrangian  L given 
by  relation  (6.27)  if^c(t)  - f ' (K(t))Jis  positive  on  (•rp)  and 
maximize  it  if(^c(t)  - f (k(t))Jis  negative  on  (»,e)  . No  variational 
principle  exists  if^c(t)  - f*(k(t))j  changes  signs.  In  that  case 
the  best  possible  statement  is  that  the  solution  of  (6.2)  will 
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rt,-  1*;^-  f 


correspond  to  a local  stationary  behavior  of  the  Lagrangian. 

Vice-vcrsa  if  we  choose  arbitrarily  p(t)  and  compute  K(x(t)) 
using  the  relation  W„  » ^*p,  then  the  choice  of  p which  will  minimize 

J\ 

the  Lagrangian  corresponds  to  the  solution  of  the  system  (6.2). 

Since  oAly  once  differentiable  functions  were  considered  in  our 
variational  arguments,  the  solutions  we  talk  about  are  weak  solutions 
in  Sobolev  space  (or  («<iB)  in  a commonly  used  symbolism.) 
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